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Abstract 

Linear cosmological perturbation theory is pivotal to a theoretical understanding of current cosmological 
experimental data provided e.g. by cosmic microwave anisotropy probes. A key issue in that theory is to extract 
the gauge invariant degrees of freedom which allow unambiguous comparison between theory and experiment. 

When one goes beyond first (linear) order, the task of writing the Einstein equations expanded to n'th order 
in terms of quantities that are gauge invariant up to terms of higher orders becomes highly non-trivial and 
cumbersome. This fact has prevented progress for instance on the issue of the stability of linear perturbation 
theory and is a subject of current debate in the literature. 

In this series of papers we circumvent these difficulties by passing to a manifestly gauge invariant framework. 
In other words, we only perturb gauge invariant, i.e. measurable quantities, rather than gauge variant ones. 
Thus, gauge invariance is preserved non perturbatively while we construct the perturbation theory for the 
equations of motion for the gauge invariant observables to all orders. 

In this first paper we develop the general framework which is based on a seminal paper due to Brown 
and Kuchar as well as the realtional formalism due to Rovelli. In the second, companion, paper we apply our 
general theory to FRW cosmologies and derive the deviations from the standard treatment in linear order. As it 
turns out, these deviations are negligible in the late universe, thus our theory is in agreement with the standard 
treatment. However, the real strength of our formalism is that it admits a straightforward and unambiguous, 
gauge invariant generalisation to higher orders. This will also allow us to settle the stability issue in a future 
publication. 



* gieskri@aei.mpg.de 
^ stefa n Onord ita .org 

"fthiemannOaei.mpg.de.tthiemannOperimeterinstitute.ca 
^owinklerOperimeterinstitute.ca 



1 



Contents 

1 Introduction 3 

2 The Brown — Kuchaf formalism |^ 

2.1 Lagrangian Analysis 11 

2.2 Hamiltonian Analysis 

2.3 The Brown - Kuchaf Mechanism for Dust 

2.3.1 Deparametrisation: General Theory 

2.3.2 Deparametrisation: Scalar Fields 

2.3.3 Deparametrisation: Dust 

2.3.4 Deparametrisation for Dust: Sign Issues 

2.4 Dust Interpretation 

3 Relational Observables and Physical Hamiltonian 

3.1 Implementing spatial diffeomorphism invariance 

3.2 Implementing invariance with respect to the Hamiltonian constraint 

3.3 Constants of the physical Motion 

4 Physical Equations of Motion 

4.1 First Order (Hamiltonian) Formulation 

4.2 Second Order (Lagrangian) Formulation 

4.3 Derivation of the Second Order Equations of Motion 

4.4 Summary of Second Order Equations of Motion 

5 Physical Hamiltonian, Boundary Term and ADM Hamiltonian 

6 Linear, Manifestly Gauge Invariant Perturbation Theory 

6.1 Derivation of the Equation of Motion for 63 

6.2 Derivation of the Equation of Motion for 6Qjk 

6.3 Summary of the Equations of Motion for 63 and 5(5jfc 

7 Comparison vi^ith Other Frameworks 

8 Conclusions and Open Questions 

A Second Class Constraints of the Brown — Kuchaf Theory 

B Comparison with Symplectic reduction 

C Effective Action and Fixed Point Equation 

D Two Routes to Second Time Derivatives of Linear Perturbations 

E Constants of the Motion of n'th Order Perturbation Theory 

F Generalisation to Other Deparametrising Matter 

G Linear Perturbation Theory: Some Calculations in More Detail 

H Gauge Invariant Versus Gauge Fixed Formalism 



2 



1 Introduction 



General relativity is our best theory for gravitational physics and, sofar, has stood the test of time and 
experiments. Its complicated, highly nonlinear equations of motion, however, mean that the calculation of 
many gravitational processes of interest has to rely on the use of approximations. One important class of 
such approximations is given by perturbation theory, where, generally speaking, one perturbs quantities of 
interest, such as the metric and matter degrees of freedom around an exact, known solution which, typically, 
displays a high degree of symmetry. 

It is well-known that perturbation techniques in general relativity pose challenges above and beyond 
those typically associated with them in other areas of physics, such as stability, convergence issues etc. 
The reason is that general relativity is a gauge theory, the gauge group being the diffeomorphism group 
Diff(M) of the spacetime manifold M. As a result, all metric and matter variables transform non-trivially 
under gauge transformations. This creates the problem of differentiating between (physical) perturbations 
of a given variable and the effect of a gauge-transformation on the latter. The obvious solution to this 
situation would be to calculate only with observables and perturb those. It has proved extremely difficult, 
however, to find observables in the full theory, with the exception of a few special situations, such as for 
asymptotically flat spacetimes. As a way out of this conundrum, one usually resorts to calculating in a 
specific gauge, carefully ensuring that all calculated quantities are gauge-independent. Alternatively, one 
tries to construct quantities that are observables up to a certain order. In the cosmological standard model 
this has been successfully done in linear order and forms an integral part of the modern lore of cosmology. 
In fact, there have been attempts to extend this even to second order and beyond, see, e.g., [H [21 [3l S] . The 
sheer complexity of those calculations, however, shows that there is a natural limit to how far that approach 
can be pushed. Furthermore, it is not clear whether it will succeed for other backgrounds, such as a black 
hole spacetime etc. 

This clearly makes the search for a more general framework for perturbation theory of observable quan- 
tities highly desirable. Another motivation comes from the prospects of developing perturbation methods 
for non-perturbative quantum gravity approaches, such as loop quantum gravity [5j. It is clear that the 
standard methods mentioned earlier will be extremely difficult, if not impossible to implement. 

This paper, the first in a series dedicated to this challenge, lays the foundations at the level of the full 
theory. Subsequent papers will deal with simplified cases of particular interest, such as perturbations around 
an FRW background. 

After this brief overview of the motivations behind our paper, let us now discuss some of these issues in 
more detail. The crucial ingredient in our undertaking is the construction of observables for the full theory. 
To that end let us first recall the counting of the true degrees of freedom of general relativity: The temporal- 
temporal as well as the temporal-spatial components of the Einstein equations do not contain temporal 
derivatives of four metric functions (known as lapse and shift). Thus, in the Lagrangean picture, these four 
sets of equations can be used, in principle, in order to eliminate the temporal-temporal and temporal-spatial 
components of the metric in terms of the spatial-spatial componentqH In addition, diffeomorphism gauge 
invariance displays four additional degrees of freedom as pure gaugcl. This is why general relativity in 
vacuum (without matter) has only two true (configuration) degrees of freedom (gravitons). 

In the canonical picture, the ten equations split into four plus six equations. The four equations are the 
afore mentioned constraints which canonically generate spacetime diffeomorphisms, that is, gauge transfor- 
mations. The other six equations are canonically generated by a canonical "Hamiltonian" which is actually 
a linear combination of these constraints, and thus also generates gauge transformations and even is con- 
strained to vanish. It is customary not to call it a Hamiltonian but rather a Hamiltonian constraint. The 

^In the Hamiltonian picture, these equations relate canonical momenta to canonical configuration coordinates. There are 
four additional (so called primary) constraints which impose that the momenta conjugate to lapse and shift vanish which leaves 
only two independent momenta. These eight constraints are of the first class type in Dirac's terminology [6]. 

■^In the Hamiltonian picture, the eight constraints canonically generate gauge transformations which displays eight out of 
ten configuration variables as pure gauge. 
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interpretation of Einstein's equations as evolution equations is therefore unconvincing. Instead, the correct 
interpretation seems to be that they actually describe the flow of unphysical degrees of freedom under gauge 
transformations. Thus we contend that their primary use is to extract the true degrees of freedom in the 
way described below. These true degrees of freedom are gauge invariant and thus have trivial evolution with 
respect to the canonical Hamiltonian (constraint). This is the famous problem of time of General Relativity 
[8]: There is no true Hamiltonian, only a Hamiltonian constraint and the observable quantities do not move 
under its flow. Nothing seems to move, everything is frozen, in obvious contradiction to our experience. 
This begs, of course, the question of what determines the time evolution of the true physical observables. 

In [9] a possible answer was proposed. Namely, it was shown that the problem of time can be resolved 
without affecting the interpretation of Einstein's equations as evolution equations by adding certain matter 
to the system. The method for doing this is based on Rovelli's relational formalism [10], which was recently 
extended considerably by Dittrich [IT], as well as on the Brown - Kuchaf mechanism [12]. This necessarily 
uses a canonical approach. Furthermore, it was shown in [9j that this in one stroke provides the true degrees 
of freedom and provides us with a true (physical) Hamiltonian which generates a non trivial evolution of the 
gauge invariant degrees of freedom. Remarkably, these evolution equations look very similar to Einstein's 
equations for the type of matter considered. The type of matter originally used in [9] was chosen somewhat 
ad hoc and guided more by mathematical convenience rather than physical arguments H. Furthermore, it 
seems desirable to find the optimal matter which would minimally affect the standard interpretation of 
Einstein's equations as evolution equations while increasing the number of true degrees of freedom by four. 
As it turns out, there is a natural candidate, which we will use for our purposes: Pressure free dust as 
introduced in the seminal paper by Brown and Kuchaf [T2j cited before. The dust particles fill time and 
space, they are present everywhere and at every instant of time. They follow geodesies with respect to the 
dynamical four metric under consideration. However, they only interact gravitationally, not with the other 
matter and not with itself. The dust serves as a dynamical reference frame solving Einstein's hole problem 
|14j . It can be used to build gauge invariant versions of all the other degrees of freedom. 

The dust supplies the physically meaningless spacetime coordinates with a dynamical field interpretation 
and thus solves the "problem of time" of General Relativity as outlined above. This is its only purpose. For 
every non - dust variable in the usual formalism without dust there is unique gauge invariant substitute 
in our theory. Once these observables, that is gauge invariant matter and geometry modes, have been 
constructed as complicated aggregates made out of the non gauge invariant matter, geometry and dust 
modes, the dust itself completely disappears from the screen. The observable matter and geometry modes 
are now no longer subject to constraints, rather, the constraints are replaced by conservation laws of a 
gauge invariant energy momentum density. This energy momentum density is the only trace that the dust 
leaves on the system, it can be arbitrarily small but must not vanish in order that the dust fulfills its role 
as a material reference frame of "clocks and rods" . The evolution equations of the observables is generated 
by a physical Hamiltonian which is simply the spatial integral of the energy density. These evolution 
equations, under proper field identifications, can be mapped exactly to the six of the Einstein equations for 
the unobservable matter and geometry modes without dust, up to modifications proportional to the energy 
momentum density. Thus again the influence of the dust can be tuned away arbitrarily and so it plays a 
perfect role as a "test observer" . It is interesting that in contrast to |I2| the dust must be a "phantom dust" , 
for the same reason that the phantom scalars apperared in [9J : If we would use usual dust as in [12] then 
the physical Hamiltonian would come out negative definite rather than positive definite. Or equivalently, 
physical time would run backwards rather than forward. Notice that general relativistic energy conditions 
for the gauge invariant energy momentum tensor are not violated because it does not contain the dust 
variables and it is the dust free and gauge invariant energy momentum tensor that the positive physical 
Hamiltonian generates. Hence, while the energy conditions for the phantom dust species are violated at 
the gauge variant level, at the gauge invariant level there is no problem because the dust has disappeared. 
Notice also that even at the gauge variant level the energy conditions for the total energy momentum tensor 

''Also, apart from cosmological settings, the consequences of the deviations of these evolution equations from Einstein's 
equations was not analysed. 
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are still satisfied if there is sufficient additional, observable matter present. 

Based on these constructions we will develop a general relativistic perturbation theory in this series of 
papers. In the current work we treat the case of a general background, in the follow-up papers we discuss 
special cases of particular interest. 

The plan of this paper is as follows: 

In section [2] we review the seminal work of Brown and Kuchaf [12]. We start from their Lagrangian 
(with opposite sign in order to get phantom dust) and then perform the Legendre transform. This leads to 
second class constraints which were not discussed in [12] and which we solve in appendix [XI After having 
solved the second class constraints the further analysis agrees with [12]. The Brown - Kuchaf mechanism 
can now be applied to the dust plus geometry plus other matter system and enables us to rewrite the four 
initial value constraints of General Relativity in an equivalent way such that these constraints are not only 
mutually Poisson commuting but also that the system deparametrises. That is, they can be solved for the 
four dust momentum densities, and the Hamiltonian densities to which they are equal no longer depend on 
the dust variables. 

In section [3] we pass to the gauge invariant observables and the physical Hamiltonian. In situations 
such as ours where the system deparametrises, the general framework of [11] drastically simplifies and 
one readily obtains the Dirac observables and the physical Hamiltonian. Due to general properties of the 
relational approach, the Poisson algebra among the observables remains simple. More precisely, for every 
gauge variant non - dust variable we obtain a gauge invariant analog and the gauge variant and gauge 
invariant observables satisfy the same Poisson algebra. This is also proved for part of the gauge invariance 
by independent methods in appendix [Bj The physical time evolution of these observables is generated by a 
unique, positive Hamiltonian. 

In section S] we derive the equations of motion generated by the physical Hamiltonian for the physical 
configuration and momentum observables. We also derive the second order in time equations of motion 
for the configuration observables. Interestingly, these equations can be seen of almost precisely the usual 
form that they have in the canonical approach [26] if one identifies lapse and shift with certain functions of 
the canonical variables. Hence we obtain a dynamical lapse and shift. The system of evolution equations 
is supplemented by four sets of conservation laws which follow from the mutual commutativity of the 
constraints. They play a role quite similar to the initial value constraints for the system without dust 
written in gauge variant variables but now the constraint functions do not vanish bur rather are constants 
of the motion. 

In section [5] we treat the case of asymptotically flat spacetimes and derive the necessary boundary terms 
to make the Hamiltonian functionally differentiable in that case. Not surprisingly, the boundary term is 
just the ADM Hamiltonian. However, while in the usual formalism the bulk term is a linear combination 
of constraints, in our formalism the bulk term does not vanish on the constraint surface, it represents the 
total dust energy. 

In appendix [Owe perform the inverse Legendre transform from the physical Hamiltonian to an action. 
This cannot be done in closed form, however, we can write the transform in the form of a fix point equation 
which can be treated iteratively. The zeroth iteration precisely becomes the Einstein - Hilbert action 
for geometry and non - dust matter. Including higher orders generates a more complicated "effective" 
action which contains arbitrarily high spatial derivatives of the gauge invariant variables but only first time 
derivatives. 

In section [6] we perturb the equations of motion about a general exact solution to first order, both in the 
first time derivative order form and the second time derivative order form. Notice that our perturbations 
are fully gauge invariant. In appendix [D] we show that one can get the second time derivative equation of 
motion for the perturbations in two equivalent ways: Perturbing the second time derivative equations of 
motion to first order or deriving the second time order equation from the perturbations to first order of the 
first time order equations. This is an important check when one derives the equations of motion for the 
perturbations on a general background and the second avenue is easier at linear order. However, the first 
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avenue is more economic at higher orders. In appendix [E] we show that the equations of motion up to n'th 
order are generated by the physical Hamiltonian expanded up to (n+l)th order. Moreover we show that 
the invariants expanded to n'th order remain constants of the motion under the (n+l)th order Hamiltonian 
up to terms of at least order n+1. This is important in order to actually derive the second time derivative 
equations of motion because we can drop otherwise complicated expressions. 

In section eight we compare our new approach to general-relativistic perturbation theory with some 
other approaches that can be found in the literature. 

Finally, in section nine we conclude and discuss the implications and open problems raised by the present 
paper. 

In appendix |F] we ask the question whether the qualitative conclusions of the present paper are generic or 
whether they are special for the dust we chose. In order to test this question we sketch the repetition of the 
analysis carried out for the phantom dust for the phantom scalar field of [9]. It seems that qualitatively not 
much changes, although the dust comes closer than the phantom scalar to reproducing Einstein's equations 
of motion. This indicates that the Brown - Kuchaf mechanism generically leads to equations of motion 
for gauge invariant observables which completely resemble the equations of motion of their gauge variant 
counter parts. 

Appendix [G] contains more details concerning some calculations in section seven. 

Appendix |H] derives the connection between our manifestly gauge invariant formalism and a correspond- 
ing gauge fixed version of it. 

Finally, our rather involved notation is listed, for the convenience of the reader, on the next page. 
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Notation 



As a rule of thumb, gauge non invariant quantities are denoted by lower case letters, gauge invariant 
quantities by capital letters. The only exceptions from this rule are the dust fields T, , p, Wj, their conju- 
gate momenta P, Pj , /, P and their associated primary constraints Zj , Z, Z^ which however disappear in the 
final picture. Partially gauge invariant quantities (with respect to spatial diffeomorphisms) carry a tilde. 
Background quantities carry a bar. Our signature convention is that of relativists, that is, mostly plus. 
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symbol 



meaning 



Gn Newton constant 

K = IGttGn gravitational coupling constant 

A scalar coupling constant 

A cosmological constant 

M spacetime manifold 

X spatial manifold 

T dust time manifold 

S dust space manifold 

v, p, .. = 0, .., 3 tensor indices on M 

a, b, c, .. = 1, 2, 3 tensor indices on X 

k, .. = 1, 2, 3 tensor indices on S 

Xf^ coordinates on M 

x"" coordinates on X 
coordinates on S 

t foliation parameter 

r dust time coordinate 

Yf' one parameter family of embeddings X ^ M 

Xt = Yt{X) leaves of the foliation 

(7^1/ metric on M 

Qah (pullback) metric on X 

(jij (pullback) metric on S 

Qij Dirac observable associated to qab 

p"-^ momentum conjugate to qab 

p^^ momentum conjugate to q^j 

P^^ momentum conjugate to Qij 

C scalar field on M 

^ scalar field on X 

^ pullback scalar field on S 

H Dirac observable associated to ^ 

TT momentum conjugate to ^ 

TT momentum conjugate to ^ 

n momentum conjugate to S 

V potential of C, i-, S 

T dust time field on X 

T dust time field on S 

dust space fields on X 

p dust energy density on M, X 

Wj dust Lagrange multiplier field on M, X 
U = ~dT + WjdS^ dust deformation covector field on M 

J = det{dS/dx) dust field spatial density on X 

P momentum conjugate to T 

P momentum conjugate to T 

Pj momentum conjugate to 

I momentum conjugate to p 

P momentum conjugate to Wj 

Zj, Z, Z^ dust primary constraints on X 

p,^ , p, Pj dust primary constraint Lagrange multipliers on X 
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difFeomorphism of X 

unit normal of spacelike hypersurface on M 

n coordinate lapse function on X 

n" coordinate shift function on X 

p momentum conjugate to n 

Pa momentum conjugate to n" 

z, Za primary constraint for lapse, shift 

ly, v"" lapse and shift primary constraint Lagrange multipliers 

0, ip, B, E MFB scalars on X, S 

Sa, Fa MFB transversal vectors on X 

Sj , Fj MFB transversal vectors on S 

hah MFB transverse tracefree tensor on X 

hjk MFB transverse tracefree tensor on S 

* linear gauge invariant completions of ^, ^ 

Va linear gauge invariant completions of Fa 

Vj linear gauge invariant completions of Fj 

c*°* total spatial difFeomorphism constraint on X 

c*°* = SjC^a^ total spatial difFeomorphism constraint on X 

c*°* total Hamiltonian constraint on X 

Ca non - dust contribution to spatial diffeomorphism constraint on X 

Cj = SjCa non - dust contribution to spatial diffeomorphism constraint on X 

Cj non - dust contribution to spatial diffeomorphism constraint on S 

Cj 7^ Cj momentum density: Dirac observable associated to Cj 

c non - dust contribution to Hamiltonian constraint on X 

c non - dust contribution to Hamiltonian constraint on S 

C ^ c Dirac observable associated to c 

h energy density on X 

h energy density on S 

H = h energy density: Dirac observable associated to h 

hj = c*°* — Pj auxiliary density on X 

€ numerical energy density on S 

Cj numerical momentum density on S 

H = d?a H physical Hamiltonian, energy 

L Lagrange density associated to H 

L = d'^a L physical Lagrangian 

Vjk velocity associated to Qjk 

T velocity associated to S 

N = C/H dynamical lapse function on S 

Nj = —Cj/H dynamical shift function on S 

= Q^^N]^ dynamical shift function on S 

V ^ g^u compatible covariant differential on M 

Da Qab compatible covariant differential on X 

Dj Qjk compatible covariant differential on S 

Dj Q jk compatible covariant differential on S 

Qjk background spatial metric 

P^^ background momentum conjugate to Qjk 

S background scalar field 

n background momentum conjugate to S 

p = ^["^ + vCE)] background scalar energy density 

p = ^[S^ — v{E)] background scalar pressure 

9 



Gjkmn = Qj(mQn)k " ^QjkQmn physical DeWitt bimetric 

^Q-iykmn ^ Qj{mQn)k _ QjkQmn inverse physical DcWitt bimetric 

Gjkmn = ^j{m^n)k ~ ^^jk^mn flat background DcWitt bimetric 

^Q-ijjkmn ^ §j{mgn)k _ jifc^mn] inverse flat background DeWitt bimetric 
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2 The Brown — Kuchar formalism 



In this section we review those elements of the Brown - Kuchar formahsm [12j that are most relevant to us. 
Furthermore, we present an explicit constraint analysis for the system where gravity is coupled to a generic 
scalar field and the Brown ~ Kuchar dust, based on a canonical analysis using the full arsenal of Dirac's 
algorithm for constrained Hamiltonian systems. 

For concreteness, we employ dust to deparametrise a system consisting of a generic scalar field ( on 
a four-dimensional hyperbolic spacetime {M,g). The corresponding action, Sgeo + •S'matter! is given by the 
Einstein - Hilbert action 

5geo = - / d^X ^Ad^\ [R^"^ + 2A] (2.1) 

where k = IGttGn, with Gn denoting Newton's constant, R^^^ is the Ricci scalar of g and A denotes the 
cosmological constant, and the scalar field action 

.Sn^atter = ^ [ d^X ^/ditM ["^^'C/^C " v{C)] (2-2) 

with A denoting a coupling constant allowing for a dimensionless C, and u is a potential term. 



2.1 Lagrangian Analysis 

In their seminal paper [12] Brown and Kuchar introduced the following dust actioiill 

Sdust = -\ f d^X V\det{g)\ p [g^^-^ U^U, + 1] . (2.3) 

Here, g denotes the four-metric on the spacetime manifold M. The dust velocity field is defined by f7 = 
— dT + WjdS^ [j e 1,2,3. The action Sdust is a functional of the fields /O, g^u, T, , wj^ The physical 
interpretation of the action will now be given in a series of steps. 
First of all, the energy momentum of the dust reads 



By the Euler-Lagrange equation for p 

55dust _ gXajj^^^ ^^^Q (2.5) 



6p 

the second term in (j2.4p vanishes on shell. Hence, U is unit timelike on shell. Comparing with the energy 
momentum tensor of a perfect fiuid with energy density p, pressure p and unit (timelike) velocity field U 

T^l = P U^U, + p {g^, + U^U,) (2.6) 

shows that the action (|2.3p gives an energy-momentum tensor for a perfect fiuid with vanishing pressure. 
For p / 0, variation with respect to Wj yields an equation equivalent to 

CuS^ = (2.7) 

where C denotes the Lie derivative. Hence, the fields are constant along the integral curves of U . Equation 
(j2.7p implies 

CuT = W'T^ = [/^[T^ - WjS^^] = -U^'U^ = +1 (2.8) 



A classical particle interpretation of this action will be given in section 12.41 
^Here, T,S-' have dimension of length, Wj is dimensionless and, thus, p has dimension length"''. The notation used here is 
suggestive: T stands for time, for space and p for dust energy density. 
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so that T defines proper time along the dust flow lines. 
Variation with respect to T results in 

9^[pv^ditMf^1 = V\^^\V^[pU^] = (2.9) 

while variation with respect to gives 

d^[py'\det{g)\U^'W,] = y'l dei{g)\V ^[pU^'Wj] = . (2.10) 

Using (j2.9p . (|2.10p reduces to (assuming /> / 0) 

VuWj = . (2.11) 

Thus, VuU^ = 0, and, as a consequence, the integral curves of U are affinely parametrised geodesies. The 
physical interpretation of the fields T, is complete: the vector field U is geodesic with proper time T, 
and each integral curve is completely determined by a constant value of . This determines a dynamical 
foliation of M, with leaves characterized by constant values of T. A given integral curve intersects each leave 
at the same value of . 

2.2 Hamiltonian Analysis 

In this section we derive the constraints that restrict the phase space of the system of a generic scalar field 
on a spacetime {M,g), extended by the Brown - Kuchaf dust. The reader not interested in the details of 
the derivation, which uses the full arsenal of Dirac's algorithm for constrained systems, may directly refer 
to the result ([2M23i]) . 

We assume (M, g) to be globally hyperbolic in order to guarantee a well posed initial value problem. As a 
consequence, M is diffeomorphic to M x ^Y, where X is a three-manifold of arbitrary topology. The spacelike 
leaves of the corresponding foliation are obtained as images of a one parameter family of embeddings 
t I— > Y^, see e.g. [26] for more details and our notation table for ranges of indices etc. The timelike unit 
normals to the leaves may be writteijfl as = [Y^^ — n'^Y^a\/n, where n, are called lapse and shift 
functions, respectively. Throughout, is assumed to be future oriented with respect to the parameter t, 
which requires n > 0. 

The three metric on X is the pull back of the spacetime metric under the embeddings, that is, Qabix, t) = 
Y^oY^Qfiu- Denoting the inverse of qah by g"^ it is not difficult to see that 

gl^U ^ _^^,^U ^ ^ab _ (2.12) 

It follows that the dust action can be written as 

5dust = -\ I dt [ d^x v'det(g) n p (-U^ + q'^'UaUb + l) (2.13) 

with = n^f/^, \Ja = Y%Uf,. 

The form ()2.13p is useful to derive the momentum fields canonically conjugate to T,S^, respectively, as 

-y/det{q) p Un 



y^d^pUnWj. (2.14) 

The second relation in (12.141) shows that the Legendre transform is singular, and we obtain the primary 
constraint (Zwangsbedingung) 

Zj := Pj + WjP = . (2.15) 

^We have written Y{t,x) = Yt{x). 



6Sd 



ust 



6Si 
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Additional primary constraints arise when we compute the momenta conjugate to p and Wj 



(2.16) 



Considering the total action S = Sgeo + 'S'matter + 5'dust , further primary constraint follow from the calculation 
of the canonical momentum fields conjugate to lapse and shift n,n", respectively, 



p := z :- 

Pa ■■= Za : 



5S_ _ 
6S_ 











(2.17) 



The primary constraints signify the fact that we cannot solve for the velocities {S\, p^t: ,t-, "^,1-, ''^"f}) 
respectively, in terms of the momenta and configuration variables. Therefore, all primary constraints must 
be included in the canonical action, together with appropriate Lagrange multipliers {/x-^, p, Pj,i^, i^"}, in 
order to reproduce the Euler - Lagrange equations. 

It is straightforward to solve for and Qab ,t- For instance. 



Tt = nTn + n'^Ta = n [-Un + WjSi] + n'^Ta = n-^=£= + S\Wj + n"^ [T, - WjS\ 

P Vdet(g) 



(2.18) 



How to eliminate the velocities of the scalar field and the three- metric is well known, e.g. [26j, and will not 
be repeated here. 

The resulting Hamiltonian constraint for the extended system, c*°* = c^'^° + c™'^**'^'^ + c'^^^^^ is explicitly 
given by 



,geo 



1 



V'det(g) 



1 

Qacqhd — -^qabqcd 



p'^Y'^ - ydit(^ + 2A^/d^ 



Ac 



.matter 



„dust 



TT 



+ \/ditM P (q'^'UaUk + 1 
ydet(g) V 



(2.19) 



with Ua = —T^a + Wj S-^^a- The spatial diffeomorphism constraints for the extended system, c*°* = cf^" + 



_matter i _dust 



+ c'^ , are explicitly given by 



,geo 



-2 qacDb p" 

^dust ^ p [^^^ _ 5,^] 

The total action in canonical form reads 

d^x (pTt + Pj S\ + 1 p,t + P Wj,t + pnt+ Pa n\ + V' qah,t + It 

primary 



(2.20) 



dt 
dt H, 



(2.21) 



with p""^ denoting the momentum field conjugate to Qab, C denoting the pullback of C to A", vr denoting its 
canonical momentum and D is the covariant differential compatible with Qab- Furthermore, the Hamiltonian 
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and spatial diffeomorphism constraints, together with the primary constraints, entered the definition of the 
primary Hamiltonian 



-f^primary — / d X /iprimary (2.22) 



via the density 

Vimary = /i^' + fl Z + flj Z^ + U Z + Za + n C*°* + C*°* . (2.23) 

Consistency requires that the the constraint surface, defined by the primary constraints p. 150 . ()2.16p and 
(j2.17p . is stable under the action of -ff primary In order to to check this, we summarise the only non- vanishing 
elementary Poisson bracket^ 



{p''\x),qM} 






= A(5(3;,y) 


{P{x),T{y)} 


= Hx,y) 


{P,ix),S'iy)} 


= Hx,y) 


{iix),piy)} 


= 6{x,y) 


{P{x),Wkiy)} 


= ^i^ix,y) 


{p{x),n{y)} 


= S{x,y) 


{Pa{x),n''{y)} 





(2.24) 

The primary constraints transform under the action of the primary Hamiltonian -ffprimary as follows 



— { -f^primary ) — C 



primary) j 

tot 



^a,t — { -f^primary ) Pa } — C^j 



Z^t — {-f^primary) 1} — „ 



Y^det(g) ^ 



Z^, = {//primary, P} = - P - n p ^/ dei{q) q''' Ua ^ + P 

^ npy^det(g) 



t — 1-'-' primary 5 

Zj,t = {i/primary, Pj + W,P} = fij P - (u^ " ^^V'det(g) ^ ^ ^^^^ _ ^^.25) 

Consistency demands that (j2.25p must vanish. Indeed, the last two equations in (j2.25p involve the Lagrange 
multipliers /i-', /ij, respectively, and can be solved for them, since the system of equations has maximal rank. 
However, the first three equations in (12.250 do not involve Lagrange multipliers. Hence, they represent 
secondary constraints. According to Dirac's algorithm, the secondary constraints in equation (j2.25p force 
us to reiterate the stability analysis, i.e. to calculate the action of //primary on the secondary constraints. A 
lengthy calculation presented in |A] shows that the secondary constraints are stable under the Hamiltonian 
flow generated by //primary In other words, no tertiary constraints arise in the stability analysis for the 
secondary constraints. However, the action of //primary on c involves the Lagrange multipliers fi^ ,fj,j,n, and 
can be solved for /i. 

The final set of constraints is given by {c*°*, c*°*, c, Zj, Z\ Z, z^, z} and it remains to classify them 
into first and second class, respectively. Obviously, 

{Z^(x),Zfc(y)} = P5^,5(x,y) 

{Z(x),c(y)} = ^^ ^(x,^) (2.26) 



Notice that n,n'^ ,Wj, p, are not Lagrange multipliers at this point, they are canonical coordinates just like the other 
fields. 
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does not vanish on the constraint surface defined by the final set of constraints, hence they are of second 
class. Next, since n appears at most linearly in the constraints, while does not appear at all, it follows 
immediately that of first class. Further, consider the linear combination of constraints 

= I P,a + P Wj^a + PTa + Pj S^^, + p n,, + £^ Pa + Ca 

= +Z p,a + Wj,a + Zj S^, + zna + Za (2.27) 

where 

Ca = cr + C'**'' (2.28) 

is the non-dust contribution to the spatial diffeomorphism constraint cjf*. Since all constraints are scalar or 
covector densities of weight one and c|^°* is the generator of spatial diffeomorphisms, it follows that c|^°* is 
first class. Finally, we consider the linear combination 

gtot ^ ^tot ^ 2. + + aZ (2.29) 

and determine the phase space functions such that c*"' has vanishing Poisson brackets with 

Zj, Z^ , Z up to terms proportional to Zj, Z^ , Z. Then, c*°* is first class, as well. See appendix [X] for details. 

In the final step we should calculate the Dirac bracket O [22] {f,g}* for phase space functions f,g. It 
differs from the Poisson bracket {f,g} by linear combinations of terms of the form {f,Zj{x)} {g,Z^{y)} 
and {/, Z(x)} {g,c{y)} (and terms with f,g interchanged). Fortunately, the Dirac bracket agrees with 
the Poisson bracket on functions f,g which only involve {T, , qab, n, n"} and their conjugate momenta 
{P, Pj, P""^ , p, Pa} on which we focus our attention in what follows. Using the Dirac bracket, the second 
class constraints can be solved strongly: 

-Pj/P 



p2 

q'^^UaUb+l] . (2.30) 





= 




w, = - 


z^ 


= 




P = 


z 


= 




/ = 


c 


= 







det(g) 

From the last equation in (j2.30p we find 



V'det(g) 



{Vq'^'UaUt + 1) ' , (2.31) 



with e = ±1. We may also partially reduce the phase space subject to (I2.30p by setting z = Za = and 
treating n, as Lagrange multipliers, since they are pure gauge. Then, we are left with two constraints 

Ca' = c, + cf^* (2.32) 

where 



^geo ^ ^matter ^2.33) 



and 



P Ta + Pj S\ . (2.34) 



dust 
^a 



Equations ()2.32I[23^ are the main result of this subsection. They represent the final constraints that 
restrict the phase space of the system consisting of a generic scalar field on (M, g), extended by dust. 
The form of the dust Hamiltonian and spatial diffeomorphism constraints jc'^'i'^* ,0^"*^*!, respectively, is of 
paramount importance for utilising dust as a deparameterising system, as we will explain in the next section. 
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2.3 The Brown — Kuchaf Mechanism for Dust 



In the previous section we have shown that the canonical formulation of a classical system, originally 
described by General Relativity and a generic scalar field theory, then extended by a specific dust model, 
results in a phase space subject to the Hamiltonian and spatial diffeomorphism constraints (|2.32t - [2TM|) . The 
primary Hamiltonian, after having solved the second class constraints, is a linear combination of those final 
first class constraints (|2. 321123^ and, thus, is constrained to vanish. This holds in general, independently 
of the matter content, and is a direct consequence of the underlying spacetime diffeomorphism invariance. 

Now, observable quantities are special phase space functions, distinguished by their invariance under 
gauge transformations. In other words, their Poisson brackets with the constraints must vanish when the 
constraints hold. In particular, they have vanishing Poisson brackets with the primary Hamiltonian .ffprimary 
on the constraint surface. This is one of the many facets of the problem of time: observable quantities do not 
move with respect to the primary Hamiltonian, because the latter generates gauge transformations rather 
than physical evolution. It follows that physical evolution must be generated by a true Hamiltonian (not 
constrained to vanish, but still gauge invariant). 

In this section we address the questions how to construct a true Hamiltonian from a given Hamiltonian 
constraint, and, how to construct observable quantities (gauge invariant phase space functions). 



2.3.1 Deparametrisation: General Theory 

The manifest gauge invariant construction of a true Hamiltonian, generating physical evolution as opposed 
to mere gauge transformations, becomes particular simple when the original system under consideration can 
be extended to a system with constraints in deparametrised form. 

Consider first a general system subject to first class constraints c/. The set of canonical pairs on phase 
space split into two sets of canonical pairs {q"',Pa) and {T^,ttj), respectively, such that the constraints can 
be solved, at least locally in phase space, for the ttj. In other words, 

ci = ^ c/=7r7 + /^/(T■^;g^pJ=0. (2.35) 

Notice that, in general, the functions hj do depend on the T"^ . The first class property guarantees that the 
c/ are mutually Poisson commuting [36]. 

A system that deparametrises allows to split the set of canonical pairs into two sets of canonical pairs 
such that (1) equation (j2.35p holds globally on phase spac^, and (2) the functions hj are independent of 
the T-^. 

Property (2) implies that the functions hj are gauge invariant. Hence, any linear combination of the 
hi that is bounded from below is a suitable candidate for a true Hamiltonian in the following sense: let 
Cr = T^cj be such a linear combination, with real coefficients in the range of T^, and consider for any 
phase space function / the expression 



Of{r) 



oo ^ 

^ n! 



,n=0 



(2.36) 

Her^, the iterated Poisson bracket is inductively defined by {cr,/}(o) = f, {cT,f}{n+i) = {ct , {cr , f} (n)} ■ 
Then, Of{T) is an observable quantity. More precisely, it is a gauge invariant extension of the phase space 
function f. Furthermore, physical time translations of 0/(t) are generated by the functions hj: 

^^ = {hj,Of{T)} (2.37) 
provided that / onlj0 depends on {q"-, pa). 



*This is not the case for Klein - Gordon fields and many other scalar field theories with a canonical action that is at least 
quadratic in the ni. 

^Notice that the substitution of the phase space independent numbers by the phase space dependent combination (r — T)^ 
is performed only after the series has been calculated. 

^"This is no restriction, since the tt/ can be expressed in terms of the (g°, pa) (using (|2.35p ). and the are pure gauge. 
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The observable quantities Of{T) can also be interpreted from the point of view of choosing a physical 
gauge. Indeed, Of{T) can be interpreted as representing the value of / in the gauge = r^. 



2.3.2 Deparametrisation: Scalar Fields 

The Brown - Kuchaf mechanism relies on the observation that free scalar fields lead to deparametrisation 
of General Relativity, as we sketch below (see [9] for a detailed discussion). 

A free scalar field contributes to the spatial diffeomorphism constraint a term of the form 

^scalar ^ ^^^^ (2.38) 

and to the Hamiltonian constraint a function of vr^ and q"''^(l),a4',b, in the absence of a potential. On the 
constraint surface, defined by the spatial diffeomorphism constraint, we have the identity 

ab scalar scalar ab 

q^'cP,,<P,b = ' " = (2.39) 

with Ca denoting the contribution to the total spatial diffeomorphism constraint that is independent of the 
free scalar field. Substitution of (j2.39p into the total Hamiltonian and spatial diffeomorphism constraints 
yields the same constraint surface and gauge flow than before. In other words, the constraints with the 
substitution (j2.39p are equivalent to the original ones. However, the new total Hamiltonian constraint does 
not depend on the free scalar field cp any more. Therefore, at least locally in phase space, we can solve the 
new total Hamiltonian constraint for the momentum field vr and write locally 

c^°\x) = tt{x) + h{x) (2.40) 

where the scalar density h of weight one is independent TT,(j) and, typically, positive definite, see [9j for 
details. 

As mentioned above, the constraint (j2.40p and h{x) are mutually Poisson commuting, which guarantees 
that the physical Hamiltonian 

H := [ d^x h{x) (2.41) 
Jx 

is observable (it has vanishing Poisson brackets with the spatial diffeomorphism constraint, because h has 
density weight one). 

This is as much as the general theory goes. There are two remaining caveats: first of all, the construction 
is only local in phase space. Secondly, the construction based on a single free scalar field requires phase 
space functions that are already invariant under spatial diffeomorphisms. Only those can be completed to 



fully gauge invariant quantitie; 



reacly : 
itiefll. 



2.3.3 Deparametrisation: Dust 

Dust described by the action (j2.3p does not entirely fit into the classification scheme given in [9J and sketched 
in the last section. It is not simply based on four free scalar fields T, , but in addition leads to second 
class constraints. However, it has a clear interpretation as a system of test observers in geodesic motion, 
and circumvents the remaining caveats mentioned at the end of the last subsection as we will see. 
Recall the final form of the Hamiltonian constraint (I2.321[232]l derived in the previous section: 

ct°t = c + eP^/l + ^fHjJTb (2.42) 

with Ua = —T^a + Wj S\a- Solving the second class constraint Zj = for Wj, we find Ua = —d^^^ /P. 
Inserting the first class spatial diffeomorphism constraint c*°* = + cj^"*^*, we arrive at the equivalent 
Hamiltonian constraint 

c'-" = c + eP^jl^^^ (2.43) 
This can be circumvented by employing e.g. three more free scalars but this would be somewhat ad hoc. 
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which is already independent of T, and Pj, but still not of the form c*°* = P + h, as required for a system 
that deparametrises. 

2.3.4 Deparametrisation for Dust: Sign Issues 

In order to bring (j2.43p into the form c*°* = P + h, we have to solve a quadratic equation. Each root 
describes only one sheet of the constraint surface, unless the sign of P is somehow fixed. As we argue below, 
this freedom will be fixed by our interpretation of the dust system as a physical reference system. 

Recall that P = —p\l det((7)f7„ and = 1, = 0. In accordance with our interpretation, 

we identify T with proper time along the dust flow lines. Thus, U is timelike and future pointing, hence 
Un < 0. It follows that sgn(P) =sgn(p), so e = 1 in dOT]) . 

In [T2] the authors assume p > 0, as it is appropriate for observable dus10. In our case, however, the 
dust serves only as a tool to deparametrise the system and is, by construction, only pure gauge. Therefore, 
we relax the restriction p > 0, when solving (|2.43|) for P: 

P^ = c^- q'^^CaCk . (2.44) 

The right hand side of (j2.44p is constrained to be non - negative, albeit it is not manifestly non - negative. 
But this causes no problem, since it is sufficient to analyse the system in an arbitrarily small neighborhood 
of the constraint surface, where — q°'^CaCi) > 0. Then, 

gtot ggn(p) h, h= \/c2 - q'^^c^cb (2.45) 

is the general solution, globally defined on (the physically interesting portion of) the full phase space. 
However, c*°* is not yet of the form required by a successful deparametrisation, because of the sign function 
which also renders the constraint non - differentiable. 

In order to utilise dust for deparametrisation, the choice P < is required. Before presenting reasons 
for this choice, we stress again that the dust itself is not observable. There are three related arguments for 
the choice of P < 0: 

1. Dynamics 

The deparametrisation mechanism supplies us with a physical Hamiltonian of the form 

H = / d^xh. (2.46) 
Jx 

In the case where dust is chosen as the clock of the system, the variation of the physical Hamiltonian 
is given by 

5H = J^d^x (j^Sc- q''' ^5ca + ^ q'"'q'"'ccCd 5qab^ . (2.47) 

For P 7^ 0, then /i 7^ (in a sufficiently small neighbourhood of the constraint surface). Hence, the 
coefficients of the variations on the right hand side of (j2.47|) are non singular. Moreover, for P 7^ 0, 
also c 7^ 0, as we see from (|2.43p . In fact, using sgn(c) = — sgn(P) (from (j2.43p ) in a neighborhood of 
the constraint surface, 

£ = -sgn(P) f7^^ (2.48) 

has absolute value no less than one. 

This would be required by the usual energy conditions if the dust would be the only observable matter. However, notice 
that only the total energy momentum is subject to the energy conditions, not the individual contributions from various matter 
species. 
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Let us now compare (|2.47|) with the differential of the primary Hamiltonian constraint in the absence 
of dust: 

Hr, 



' primary 



/ d^x {nc + n''ca) (2.49) 



which is given by (lapse and shift functions are considered as Lagrange multipliers, i.e. are phase space 
independent) 

di^prirnary= / d^X {ndc + n'' dCa) ■ (2.50) 

Jx 

Comparison between ()2.47p and (j2.50p reveals that the differentials coincide, up to the additional 
term proportional to 6qab, provided we identify n := c/h as dynamical lapse and := —q°'^ci)/h as 
dynamical shift. This is promising in our aim to derive physical equations of motions for observable 
quantities which nevertheless come close to the usual Einstein equations for gauge dependent quantities. 
However, in the standard framework the lapse function is always positive, guaranteeing that the normal 
vector field is future oriented. This fact is correctly reflected in our framework only if P < 0. 

2. Kinematics 

The identification n = c/h and = —q'^^Cb/h can also be motivated as follows: 

Consider a spacetime diffeomorphism defined by i-^ {'t^'^'') •= (^(^)) 'S'"' (^)) =■ ^''(^) and let 
(r, cj-') — > Z^(t, cj) be its inverse. We can define a dynamical foliation of M by T{X) = r =const. 
hyper surf aces. The leaves Sr of that foliation are the images of S (which is the range of the 5-^) under 
the map Z at constant r. Using the identity 

5^ = Z^;, T, + Z^. S{, (2.51) 

and = 1 , S;'^ = 0, we find C/'^ = Z^r- Thus, as expected, the foliation is generated by the 

vector field U = d/dr, which is unit timelike. 

It is useful to decompose the deformation vector field U with respect to the arbitrary coordinate 
foliation that we used before: 

U'' = g^^-U, = -nf^Un + q'^' Ub . (2.52) 



From (|2T3]) and (lOTD with e = 1 we find C/„ = -^/l + (fHJaUb. Next, 

„dust „ 

Ua = = i ■ (2-53) 

On the other hand n = c/h = sgn(P)C/„ and = —Ca/h = — sgn(P) Ua- Therefore, (|2.52p can be 
written 

U^" = -sgn(P) (n + X'^ n") . (2.54) 

Hence, the sign for which n is positive yields the correct decomposition of the deformation vector field 
U in terms of lapse and shift. This calculation reveals also the geometrical origin of the identification 

n = c/h and = —q°'^Cb/h. 

As a side remark: the identity —n^ + q'^^naUi, = — 1 is an immediate consequence of the normalisation 
of the deformation vector field, g^uU^V^ = —1. That is, the deformation vector field is timelike, future 
oriented and normalised, but not normal to the leaves of the foliation that it defines. 

3. Stability and Hat spacetime limit 

Of course, we could choose P > and use —h instead of h in order to obtain equations of motion. 
However, in that case the physical Hamiltonian would be unbounded from below, leading to an unstable 
theory. Alternatively, we could stick to +h for the equations of motion, but then the r evolution would 
run backwards. 
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Moreover, since c*°* = c + P\/l + q"-''UaUb = on the constraint surface, we would have c < for 
P > 0. Since c = c^'^" + c™'^**'^'' and c™'***'"' > 0, this would enforce c^^° < 0. Hence, flat space would 
not be a solution. 

As a side remark: for Ca/h ^ 1 and P < 0, /i ~ c, while h ~ — c for P > 0. Thus, the physical 
Hamiltonian density, with respect to dust as a physical reference system, approximates the standard 
model Hamiltonian density c™'***'^'' only for P < 0. 

We emphasise again that the dust used for deperametrisation is not observable, and should not be confused 
with observable matter. It solely provides a dynamical reference frame. 

2.4 Dust Interpretation 

In this section we derive a physical interpretation of the Brown - Kuchaf action based on the geodesic 
motion of otherwise free particles |12j . 

Consider first the action for a single relativistic particle with mass m on a background g: 

Sm = -m [ ds J-g^^ . (2.55) 

The momentum conjugate to the configuration variable is given by 

^^'^ ^-gp.xpx- 

rendering the Legendre transformation singular. This is a consequence of the reparametrisation invariance 
of the action (12.55p . Hence, the system exhibits no physical Hamiltonian, but instead a primary Hamiltonian 
constraint enforcing the mass shell condition: 

C = ^(m2+5^'^P,P.) . (2.57) 

Let us proceed to the canonical formulation. In terms of the embeddings X = Yt{x\ the particle 
trajectory reads Xi^s) = y^(^)(3;(s)), so that 

X{s) = i{s) Y^t + x%s) Ya (2.58) 

where the overdot refers to differentiation with respect to the trajectory parameter s. The momenta are 
then given by 

Pa = Y'iP^ = ^ J^_ q^Jin^ + x^ 
^-gp^XPX- 

Pt = Y^P^ = (i gtt + qab (2.59) 

^-gp^XPX- 

where gu = —n^ + qatn'^nf'- We can only eliminate the spatial velocities x"". To do this set A = gtax"" = 
qatn^x"" and B = qahX°'x^ ■ Then, 

w = g^uXf'X'' = guf + 2Ai + B. (2.60) 

On the other hand, 

- ^(f'VaVh = B + 2Ai + qabn'^nH'' . (2.61) 
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Substituting w from (|2.6Up and collecting coefficients of A^B^t yields 



1 + 



^ + ^^T^- (2-62) 



_j_ Q PaVb 



Now we can solve the first equation in (j2.59p for x"". 



x^ = t {-n^±\ll+'-^] . (2.63) 



Inserting this into the second equation in (12.590 leads to a constraint of the form C = ps + h: 



C = ps- rfpa ±n^/m? + q'^^PaPb (2.64) 

while the canonical Hamiltonian is obtained from the Lagrangian in (j2.55p as 

^canon = P^.X^' -L = iC. (2.65) 

Since the constraint (I2.64p is in deparametrised form, the phase space can easily be reduced, leading to the 
reduced action 

5'reduced = ds {paX" - h) . (2.66) 



We extend this phase space by adding a canonical pair (r, m) and consider the extended action 



S, 



extended 



j ds {m-f + pax" - h) (2.67) 



where the particle mass m is now considered as a dynamical variable. The equations of motion for m, r give 
m = and f = tyj—w. Thus, the mass is constant and r is the proper time (in the gauge s = t). 

We generalize our results now to the case of many particles. More precisely, let 5 be a label set and 
consider a relativistic particle for each label a £ S. This amounts to provide each variable appearing in the 
extended action with a corresponding label, i.e. x*^, p'^, Tu, m" ^ and the total action for those particles is 
then just the sum over the corresponding actions 'Sg^tended 

5'extended = '^'extended • (2.68) 

Next we consider the limit in which S becomes a three - manifold, with the labels a becoming coordinates 
on this manifold. In this limit, we introduce the following fields: 





fia) 




P[ 


» dV 


= nf 








Pa 


{<j)d^(J 


= Pa{Xa) 




n(o-) 


= n{xa) 




n"(fT) 


= W'ix^) 




qab{<y) 


= QabiXa] 



(2.69) 



21 



Then, in the specified hmit, the extended action (j2.68p becomes 



extended 



dt I dV (fP + + n^Pa Tn\/P^+ q'^'^PaPb 



(2.70) 



Finahy, we perform a canonical transformation: instead of the fields S"'{a) with values in X, we would 
like to consider the inverse fields S^{x) with values in S, that is S^{S{a)) = , S"'{S{x)) = x". This is at 
the same time a diffeomorphism and we can transform the other fields as well. For instance (T is a scalar 
and P is a scalar density), 



T{S{x)) 
P 



d-'a 6 (x,S{a)) det{dS/da) T{a 



-(Six)) 



det{dS/da)\ Js 

d^a aU{x,S{a)) det{dS/da 



d^a5{x,Sia)) P{a 



T{x) 
P{x) 
S^{x) 

Calculating the time derivatives and performing integrations by parts, we find 

dVrP= / d^x (TP- S^S^PTr, 



IS J X 

with denoting the inverse of the matrix S'^^a- Using 



S{X) = (7 



and defining Pj{x) implicitly through 



Pa 



PTa + PjS^a 



\det{dS/dx) 



(2.71) 



(2.72) 



(2.73) 



(2.74) 



S{x)=a 

we find that S'extended precisely turns into the dust action on X with the second class constraints eliminatecf^. 



3 Relational Observables and Physical Hamiltonian 

In this section we present an explicit prescription for constructing gauge invariant completions of arbitrary 
phase space functions. The construction is non - perturbative and technically involved, but the physical 
picture behind it will become crystal clear. Furthermore, the formal expressions are only required to establish 
certain properties of the construction, but are not required for the calculation of physical properties. This 
is a great strength of the relational formalism. 

Let us summarise the situation. After having solved the second class constraints and having identified 
lapse and shift fields as Lagrange multipliers, we are left with the following canonical pairs 



subject to the following first class constraints 



(3.1) 



tot 



„tot 



I „dust 
<-a <-a 



c + d 



dust 



dust 



„dust 



PTa + PjS^^, 



p2 _j_ ^afegdustgdust 



(3.2) 



The metric field has to be pulled back by the dynamical spatial diffeomorphism, as well. For details, see next section. 
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where Cq, c are independent of the dust variables {T, P, S^,Pj^. We already used P < 0. 

As explained in 12.3^ we aim at deparametrisation of the theory and therefore solve (j3.2p for the dust 
momenta, leading to the equivalent form of the constraints 



ct°t = P + /l , h= x/c2 - q'^'^CaCb 

cf = Pj + hj , hj = i-hTa + Ca) (3.3) 

with SjS'^^a ~ ^j^ '^j ^''b ~ hence Sj is the inverse of 5-',a (assuming, as before, that S : X ^ 
5 is a diffeomorphism) . These constraints are mutually Poisson commutinj^. However, only c*°* is in 
deparametrised form (i.e. h is independent of T,S^), but c^°* is not. In particular, we can only conclude 
that the h{x) are mutually Poisson commuting. Still, this will be enough for our purpose^. 

Following the works [QlllOtITT]. we describe the construction of fully gauge invariant completions of phase 
space functions. Consider the smeared constraint 

Kf,= [ d^x [fiix)c"'\x)+p\x)cf\x)] (3.4) 

where /3(x), (3^{x) are phase space independent smearing functions in the range of T{x), S^{x). Under a 
gauge transformation generated by this constraint, an arbitrary phase space function / is mapped to: 

oo 
n=0 

The fully gauge invariant completion of / is given by 

Oy[r,a] = U(/)] . (3.6) 

Here, the functions t{x), a^{x) are also in the range of T{x), S^{x), respectively!^. It is important to 
first calculate the Poisson brackets appearing in (|3.5p with the phase space independent functions (3,(3^, 
and afterwards to replace them with the phase space dependent functions t — T, — , respectively. 
This connection can be established based on the gauge transformation properties of T,S^: otpiT) = T + 
/?, ap{S^) = Hence, 

Of[r,a] = apif) ,^(t)=. ■ (3.7) 

Indeed, ()3.7p motivates the following interpretation: 0/[r, cr] is the gauge invariant completion of /, which 
in the gauge T = r, = takes the value /. This is not the only interpretation we entertain a different 
one below. 

For the purpose of this paper it suffices to consider the infinite series appearing in the gauge invariant 
completions as expressions useful for formal manipulations. There is no need to actually calculate these 
series for any physical problem. 

Further important properties [36\ [20] of the completion are: 

{Of[T, a],Of, [r, a]} = {0/[r, a],Of, [r, a]}* = 0{/,/'}. [r, a] (3.8) 
Of+f [r, a] = Of [r, a] + Of [r, a] , Of.f [r, a] = Of [r, a] • Of [r, a] . (3.9) 



^''One can either prove this by direct calculation, or one uses the following simple argument: the Poisson bracket between the 
constraints must be proportional to a linear combination of constraints, because the constraint algebra is first class. Since the 
constraints are linear in the dust momenta, the result of the Poisson bracket calculation no longer depends on them. Therefore, 
the coefficients of proportionality must vanish. 

In what follows we will drop the tilde in noting the constraints for notational simplicity. 

^^We denote the functional dependence of (|3.6p on the functions t{x), a-'{x) by square brackets. Below we show that it is 
sufficient to choose those functions to be constant and replace the square brackets by round ones for notational convenience. 
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Here, {•,.}* is the Dirac bracket^ [B] associated with the constraints and the gauge fixing functions T,SK 
Relations (j3.8p and (|3.9p show that the map / ^ Of[T,a\ is a Poisson homomorphism of the algebra 
of functions on phase space with pointwise multiplication, equipped with the Dirac bracket!^ as Poisson 
structure. 

In particular, for a general functional / = f[qah{x),P°-^{x),^{x),T:{x),T{x),P{x),S^ {x),Pj{x)] the fol- 
lowing useful identity holds: 

Of = /[<^ga6(2^)'^J'"Ha:)'^5(2:)'^'^(a;)'^T(x);Op(x))05j(2.),Op^,(a,)](r,Cj) . (3-11) 

This has important consequences: (jS.lip ensures that it suffices to know the completions of the elementary 
phase space variables. In fact, we are only interested in those functions that are independent of the dust 
variables {T, S'-' , P, Pj}. The reason for this is that, first of all, P, Pj are expressible in terms of all other 
variables on the constraint surface. Alternatively, since the constraints are mutually Poisson commuting, 
we have 

Op[x) = Octot(^) + 0/j(a,) = C^°^{x) + 0/j(2,) , Op^(x) = C'ctot(^) + 0^^,(2;) = C*°*(x) + Oh^(x) . (3.12) 

Hence, these functions are known, once we know the completion of the remaining variables. Secondly, 

Ot{x) [t, 0-] = t{x) , Osj(^) [r, fj] = (x) (3.13) 

are phase space independent. Thus, the only interesting variables to consider are {gab, P'^'', ^, tt}. 
In what follows we consider only dust - independent functions /. For those (j3.8p simplifies to 



{Of[r,alOf,[r,a]] = 0{fj,^[T,a]. (3.14) 

Equations (13. 9p and (13. Sp imply that / 0/[r, cr] is a Poisson automorpfiism of the Poisson subalgebra of 
functions that do not depend on the dust variables with the ordinary Poisson bracket as Poisson structure. 
This will be absolutely crucial for all what follows. 
Further useful properties of the completion are: 

0f[T,a] = 0'^^l^^^[T] (3.15) 



where (recah ([331), ([33]) and l[3M ) 



0^0 



0f[T] = K(/)],^.-T. (3.16) 

This follows from the fact that the constraints are mutually Poisson commuting and {c*'°*(x), S'-' (y)} = 0. 
The important consequence of (13.15P is that we can accomplish full gauge invariance in two stages: we 
establish first invariance under the action of the spatial diffeomorphism constraint, and afterwards achieve 
invariance with respect to gauge transformations generated by the Hamiltonian constraint. This holds even 
under more general circumstances |11) . i.e. when the constraints can not be deparametrised. 



For completeness, we note the definition of the Dirac bracket: 

{/,/'}* = {/,/'} - / ^ [{f,c\T{x)}{f,S''{x)} - {f,^°\x)}{f,S>'{x)}] (3.10) 



M = 



where 4°* = c*°', S° := T. The Dirac bracket is antisymmetric and {/, cj;°'(a;)}* = {f,S^{x)Y = everywhere. This follows 
from the fact that c^^^{x) and S^{x) are mutually Poisson commuting, and that {cj,°*(a::), = S[x,y)5'^. 

^*The Dirac bracket coincides with the Poisson bracket on gauge invariant functions. It is degenerate, since it annihilates 
coinstraints and gauge Exing functions. Hence, it defines only a Poisson structure, but not a symplectic structure on the full 
phase space. 
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3.1 Implementing spatial diffeomorphism invariance 

Keeping the physical interpretation of the completion in mind, the map / ^ oj^^[cr] can be worked out 
explicitly. In the first stage of the construction, the corresponding smeared constraint reads 

Kp = J d^x 13^ (x) cf\x) . (3.17) 

Given a phase space function /, its completion oj^^ with respect to gauge transformations generated by Kf^ 
becomes 



oo ^ 



n=0 



n=l •' 



X 



X {4f(xi),{c*f(x2),...,{4f(x„), /}...}} . 

Let us begin with / = ^(x). We claim that 

{if^,e(x)}(„) = [/3^^../3^"^,,...^;,„-e] (x) 
where Vj is the vector field defined by 

In fact the vectors are mutually commuting. 



[Wj,fA;] = SjSi,„db- j ^ k 

= -S'^S\s'Sld,-j^k 



= -sp\s\,,sid,-j^k 
= spls\,sid,-j^k 

= Sl,Sldb-j^k 

= stsl,d,-j^k 

= 

To prove (j3.19p by induction over n we need 

{Kp,S'^{x)} = -St{x){Kp,S%x)}S]{x) = -St{x)P%x)S]{x) = -[v, ■ f3^]St 
For n = 1 we have 

which coincides with (jS.lOp . Suppose that (|3.19p is correct up to n, then 



(3^^ . . . 15^^ 
. . . 



1=1 



(3.18) 
(3.19) 
(3.20) 



(3.21) 

(3.22) 
(3.23) 



■ ■ ■ ^jr. ■ - C-^Vji--- ^ji-i ['^jiP^"^'] ^h+i^ji+i ■■■yjn-i 



1=1 

n 



P^^ . 


. . p^" 


■ 






p^^ . 


..pi- 


Vj, . . 






p^^ . 




■ 







1 = 1 



(3.24) 
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where we used commutativity of the Vj and the Leibniz rule. 
It follows that 



oo 

ofl) M = + Y.-1 [^'' (^) - (^)] • • • [^'" (^) - (^)] • • • • ^(^) • (3-25) 

n=l 

Using ■ = 6j and commutativity of the Vj, we find with := — that 

oo p 
n=l L \ /• 

oo ^ 

= ^'fc • e + [I'fc • • c + X] —\^^^ ■ ■ ■ [ W ■ P^] vjvj, ...vj„-^ + vkVj, . . . ■ e] 

n=l 

oo ^ 

= Er^^''---^'" (3.26) 

n=0 



The interpretation of p.25p becomes clear for the choice a^{x) = =const., for which (j3.26p vanishes 



identically. In other words, the completion Of^^.Jf'"] does not depend on x at all. Hence, for this choice 



of , we are free to choose x in 0^|],)[o"] in order to simplify ()3.25p . Since ()3.25p is a power expansion 

in [a^{x) — S^) {x\ and is a diffeomorphism, we choose x = x^, with x^j being the unique solution of 
Si{x) = a^. Ther " 



Jl9| 



o 



\W'^=^M = U^)]s^i.)=a^ ■ (3-27) 



The completion O^l^.) (f) of ^ {x) has also a simple integral representation: 



0fl^{a) = / d\ \dei{dS{x)/dx)\ 5{S{x),(t) i{x) . (3.28) 

The significance of choosing =const is the following: Clearly, the choice a^{x) = const, is not in the 
range of S^{x), which is supposed to be a diffeomorphism. Thus, the interpretation of 0^^\a) as the value 
of / in the gauge = is obsolete. However, given a function cr^{x), instead of solving S^{x) = (x) 
for the values of the function for all x, we could solve it for x, while keeping the function arbitrary. 
This is the appropriate interpretation of 0^p{(T). This is possible because Oj^^[o"] is (at least formally) 
gauge invariant, whether or not = a-' is a good choice of gauge. It is fully sufficient to do this because, 
as shown in [12] and as we will show in appendix [Bl the partially reduced phase space (with respect to 

the spatial diffeomorphism constraint) is completely determined by the 0^p{(T), hence the 0^p[a] must be 
hugely redundant. 

We can now compute the spatially diffeomorphism invariant extensions for the remaining phase space 
variables without any additional effort, by switching first to variables which are spatial scalars on X, using 
J := det{dS/dx), which we assume to be positive (orientation preserving diffeomorphism): 

(e,7r/J) , (r,P/J) , (g,, ^ S^Sl ,pi^ = S^.S^p'^'/j) . (3.29) 

The image of these quantities, evaluated at x, under the completion O^p {a) simply consists in replacing x by 
Xct, where x^ solves S^{x) = , just as in (|3.27p . The scalars (j3.29p on X are the pull backs of the original 
tensor (densities) under the diffeomorphism a Xa- evaluated at a. Thus, they are tensor (densities) of the 
same type, but live now on the dust space manifold S. 



'We switched to the notation O'^P (a) to indicate the choice a-' (x) — a-' =const. 
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This statement sounds contradictory because of the following subtlety: We have e.g. the three quantities 
P(x), P{x) = P{x) / J{x), P{cr) = P{xa)- On P{x) is a scalar density while P{x) is a scalar. Pulling 
back P{x) to S = S{X) by the diffeomorphism a i-^ S'^^(cr) results in P{cr). But pulling back P(x) back to 
S results in the same quantity P{(t). Since a diffeomorphism does not change the density weight, we would 
get the contradiction that P{(t) has both density weights zero and one on S. The resolution of the puzzle 
is that what determines the density weight of P{x) on X is its transformation behaviour under canonical 
transformations generated by the total spatial diffeomorphism constraint c*°* = c^"'^* + where cj^"*^* , are 
the dust and non dust contributions respectively. After the reduction of c*°*, what determines the density 
weight of P{a) on S is its transformation behaviour under {[ca + PT^a]Sj / J){xa) = Cj{a) + P{a)Tj{a) and 
this shows that P{(t) has density weight onJ^ 



We will denote the images under / i— > O^^^ (a) by 



(3.30) 



In appendix iBl we show that the quantities ()3.30p can be also obtained through symplectic reduction which 
is an alternative method to show that the pairs in (|3.30p are conjugate and as it was done in |12j . 



3.2 Implementing invariance with respect to the Hamiltonian constraint 

Having completed the elementary phase space variables with respect to the spatial diffeomorphism constraint, 
it remains to render those variables invariant under the action of the Hamiltonian constraint. This amounts 
to calculate the image of those variables under the map / Oj [a] , for any / in (j3.30p . For / independent 
of T, P, the completion of / with respect to the Hamiltonian constraint is given by 

oo „ 

= E - {Mt],/}(„) , h[T] = / d'x (T(x)-r(x)) h{x) . 



(3.31) 



Only if we choose r(x) = r = const. ()3.3ip is invariant under diffeomorphisms. Hence we choose t{x) = t = 
const, which allows to rewrite (j3.3ip entirely in terms of the variables p.30p . As a reminder of this choice, 
we denote the completion by 0^p{T). In this case (j3.3ip can be written as 



with h{o-) denoting the image of h{x) under the replacement!^ of \^^{x),TT{x),qab{x),p"'''{x)^ by 
|f (cj), vf (cr), gjfc(cj),p''^(c7)|, respectively. Explicitly, denoting 



(3.32) 



c(a) 



c(x) 



S{x)=a- 



Cjjx) 
J{x) 



S{x)=a 



(3.33) 



where, as before, Cj{x) = S'-{x) Ca{x), we find 



Ha) 



(3.34) 



In order to avoid confusion of the reader we mention that any quantity f on X which has positive density weight is mapped 
to zero under / y—> Of\a). Let us again consider the example f = P. We have P{a) = P{xcr) det{dS~^ (a) /da) which is 

P(cr) = det{da/dx) P{a) = since 



J{x) P(x) 



= o 



(1) 



perfectly finite. However by the Poisson automorphism formula Op^^^ 
a = const. 

The proof of this statement is based on the fact that the replacement corresponds to a diffeomorphism and that h{r) is 
the integral of a scalar density of weight on, for r =const. 
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It is easy to see that 

±0f\T) = {H,0f\r)} (3.35) 

with 

H := / d^a h{a) (3.36) 
Js 

is the physical Hamiltonian (not Hamiltonian density) of the deparametrised system. 

We denote the fully gauge invariant completions of the Hamiltonian constraint, the spatial diffeomor- 
phism constraint^^ and the physical Hamilton density, respectively, as 

Cir,a) ^ 0|i)(r) C,(r, a) ^ og|^)(r) , 

Hir,a) ^ 0f(i)(r). (3.37) 

It is worth emphasising again that H{t, a) is the physical energy density associated to the physical Hamil- 
tonian when the dust fields are considered as clocks of the system. The fully gauge invariant completions of 
the phase space variables for matter and gravity are denoted by 

Qvir,^) - 0g(^)(r) P''^-(r,a) = 0g(^)(r). (3.38) 

The matter scalar field H(r, o") and its conjugate momentum n(T, u) are observable quantities since gauge 
invariant. The same applies to the three-metric Qij{T,a) and its canonical momentum field P^^{T,a). 
Moreover, the completion is non - perturbative, i.e. full non-Abelian gauge invariance has been accomplished. 

3.3 Constants of the physical Motion 

In the previous section we successfully constructed fully gauge invariant quantities for a specific deparmetris- 
ing system. In some sense, the construction frees the true degrees of freedom from the constraints, replacing 
them by conservation laws which govern the physical motion of observable quantities. Indeed, we have the 
following first integrals of physical motion {conservation laws): 

^Cj (r, a) = , ^H{T,a) = 0. (3.39) 
dr dr 

These equations express invariance under the physical evolution generated by H, as opposed to gauge 
invariance. The functions Cj, H, representing physical three - momentum and energy, are already gauge 
invariant. 

We proceed the proof of (j3.39p . Recall that the original constraints c*'°*(x) , c*°*(a;) are mutually Poisson 
commuting. Using (j3.3p . this means in particular, 

{c'°\x),c'°\y)} = {Fix) + hix),P{y) + hiy)} = {h{x), hiy)} = (3.40) 

where we used that the P{x) are mutually Poisson commuting and that h{x) is independent of the dust 
variables. Next, consider the smeared spatial diffeomorphism generator 



c(n) = / d'^x u'^ix) Caix) . (3.41) 
Jx 

The smeared constraint acts on h{y) as it should, 

{c{u),h{y)} = [n-h]Jy) (3.42) 



Explicit expression for the constraints in terms of the fuUy gauge invariant phase space variables are given in the next 
section, see (I3.56p . 
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or, after functional differentiation with respect to the smearing functions 

{ca{x),h{y)} = dya {S{x,y)h{y)) . 



(3.43) 



This follows from the properties of Cq, generating spatial diffeomorphims on the matter and gravity variables, 
and /i, being a scalar density of weight one and only depending on the non - dust variables. Furthermore, the 
spatial diffeomorphsims form an algebra with {c(u), c(n')} = c([u', u]). From this follows again by functional 
differentiation 

{ca{x),Cb{y)} = [dyb5{x,y)\ Ca{y) - [dx^5{y,x)] Cb{x) . (3.44) 
Let us investigate the implications of (|3. 40113^^1]) for 



h{x) 

Cgjx) S^jx) 
J{x) 



S{x)=a JX 



A^xb{S{x),o) h{x) 



S{x)=a 



X 



d^X(5(5(x),c7) S'^AxYaix) . 



First of all. 



[hia),hia')] = d^x d?y5{S{x),a) 5{S{y),a') {/i(x),%)} = 



(3.45) 



(3.46) 



where we used that the S^{x) are mutually commuting, as well as with the h{y). Second, denoting the 
puUback of the smeared diffeomorphism generator with c{u) for some smearing functions u^{(t), we have 



X 



d'x / d^y5{S{x),a)5{S{y),a') S;;{x) {caix),h{y)} 



X 



d'au^ia) / d'x / d^y6{S{x),a) 6{S{y),a') S^{x)dya{d{x,y)h{y)) 
Jx Jx 



d^au^ia) / d^x6{S{x),a) [d^a6 {S{x),a')] S'}{x)h{x) 
Jx 

d'au^ia) l'^d'x6{S{x),a) {a,a')]-^^^^,^^ S';,{x)S]{x) h{x) 

h{x) 



X 



d^au^ia) / dVi(5(ai,c7) d,5{ai,a') 
Js L 

I dVn^(ai) [<5(^Ti,c7')]^^, h{ai) 
J s 

7i u^{ai) h{ai) 5 ((Ti,cr') 
u^{a') h{a') 



J{x) 



S(x)=(n 



The last implication follows from 



c{u) = c{us) , u%{x) = S^{x)u^{S{x)) 



(3.47) 



(3.48) 



where the vector fields us are phase space dependent (they depend on S) and using the fact that the (x) 
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and Ca{y) are mutually Poisson commuting. Then, 

{c{u),c{u)} =c{[u's,us]) 



X 



d^x S^x) 



X 



u'^iSix)) [St{x)S[,{x) n5(5(x)) - Snx)Sl^ix)Sl{x)u'{S{x)) 



u^iSix)) ( Stix)u'KSix)) - S^ix)SUx)SUx)u''iS{x)) 



Ik, 



Ca{x) 



d-^x 



X 



u'^{S{x))u^^^{S{x)) - u^S{x))u[^{S{x)) Ca{x)Si{x) 



Ik, 



c{[u,u\) . 



(3.49) 



Hence, equations (jOOHOil) are exactly reproduced by (fOG]) . (fH^iTl) and ([3^ . 

We can now easily finish the proof of (j3.39p . In (j3.32p we introduced /i(t). From (j3.46p follows that 



This implies in particular that 



{h{T)M^)]=^ 



(3.50) 



(3.51) 



is already an observable quantit}o. Hence, from the definition of H and (j3.46p we find {H,h{a)} = 0. 
Furthermore, 



{H,C,(r,a)} = {0^^\r),0^>{r)} = O^'Ar) = . 



(2) 



1(2) 



{H,e,M}^ 



(3.52) 



Alternatively, a more direct way to understand this result is to make use of the series representation 
(13321) and of 



h{T)=TH-h[f] h[f] 



dV f{a) h{a) 



(3.53) 



Since the Hamiltonian vector fields Xi, X2 of H and h\T], respectively, are commuting, we may write for 
([332]) 



Cj{a,T) = exp(TXi - X2) ■ Cj{a) = exp(-X2) • [exp(rXi) • Cj{a)] 



n=0 



(3.54) 



which is clearly r - independent. 

We end this section by giving an explicit expressions for the physical Hamiltonian in terms of purely 
gauge invariant quantities: 

H{a) = ^C{T,aY-Qi\T,a) Cj{a) Ck{a) (3.55) 

Note that C, Q^'' are not independent of the physical time r. Of course, C{T,a), Cj{a) are obtained from 
c((t), Cj{a) simply by replacing everywhere the functional dependence on ^i{cr),TT{a),qjkicr),p^''{a)^ by 



■^^Note, however, although H{a) — h{a), this is not true for the corresponding spatial diffeomorphism constraints, Cj{a) 7^ 
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that on \^E{T,a),Il{T,a),Qjk{T,a), P^^{T,a)^. In greater detail, 



2Qjk{DkP''')+U{DjE) {T,a) 



Vdet(Q) V ^ J 



1 

+ 2A 



n2 



+ 



ydit(Q) (q^'^ (DjE) (DkE) + v{E) 



Vd^t(Q) 

= Cgeo {T,Cr) + Cmatter (t, CT ) 

with Dj denoting the covariant differential compatible with Qjk- 



(3.56) 



4 Physical Equations of Motion 

In this sectioij^ we derive the physical evolution of the gauge invariant functions {E, 11, Qjk, P-^'^}, 
generated by the true Hamiltonian H, in the first order (Hamilton) and second order (Lagrange) formulation. 
In other words, we study the true evolution of matter degrees of freedom and gravity with respect to the 
physical reference system (dust). 

4.1 First Order (Hamiltonian) Formulation 

For a generic observable F, we denot^^ its r-derivative simply by an overdot, F. Then, 



{H,F} = d'a {Hia),F} 
1 



(4.1) 



dV 



Cia) {Cia),F} - Q^\<j)Ck{a) {Cj{a),F} + ^g*™((T)C^(a)Q^"(a)C„(a) {Qi,(a), F} 



1 



,5 H{a) ^-----^'-^ -V , 2 

Let us introduce dynamical shift and dynamical lapse fields by 

-C,/H 



N 



C/H = Jl + Qi^NjN, 



k ■ 



(4.2) 



Notice that Nj is a constant of the physical motion, but neither are nor A^-^ = Q^^N^. Then (j4.ip can be 
rewritten in the familiar looking form 

dV (N{a){C{a),F}+N^{a){C,{a),F} + ^H{a)N\a)N\a){Q,,{a),F}] . (4.3) 



The first two terms in ()4.3p are exactly the same as those in the gauge variant derivation of the equation of 
motion, derived with respect the primary Hamiltonian 



Hprimary(iV,iV) = d^cj {N {a)C (a) + {a)Cj (a)) 



(4.4) 



For the purposes of this section we assume that X and, equivalently, S have no boundary. In order to allow for more 
general topologies, we consider boundary terms in the next section. The calculations of the present section are not affected by 
the presence of such a boundary term, because it only cancels the boundary term that would appear in the calculation of this 
section. 

Furthermore, for notational ease we drop the dependence on (r, a) when no confusions can arise. 
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Here, A'^, are viewed as phase space independent functions. The third term in (j4.3p . on the other hand, 
is a genuine correction to the gauge variant formahsm. However, it enters only in the physical evolution 
equation of P^^. Hence, 

^ N 



n 



Qjk 



N 



2N 



7dit(Q) 



N 



Qmn 

VdetQ 



2pjmpkn _ pjkpmn^ ^ I^Qjk ^ _ ^ditQ Q^'' (2A + ^ + v{E))^ 



-'^HQ^^^Q'^Nr^Nn + iCj^Py" 



(4.5) 



with Cj^ denoting the Lie derivativj^ with respect to the vector field N with components = Q^^N^, 



and we have defined the DeWitt metric on symmetric tensors as 



Gjkmn — 2 i.Q jmQnk ~l" Q jnQmk 



QjkQmn) 



which has the inverse 



.lijkmn 



QjniQnk ^ QjnQ 



anp.mk 



2QJkQmn 



(4.6) 



(4.7) 



that is Gjkmn [G ij"™^"^ = 5^5^^- The Ricci tensor of Q is denoted by Rjk[Q] and G = Cgeo+Cmatter denotes 
the split of the Hamiltonian constraint (with the dust reference system excluded) into gravitational^ and 
matter contribution, as shown explicitly in (j3.56p . 

It is already evident that the dust model we utilised as a physical reference system has the great advan- 
tage that, remarkably, equations (j4.5p are almost exactly of the same form as the corresponding equations 
in the gauge variant formalism, the only difference being the last term on the right hand side of the physical 
evolution equation for P^^. In other words, introducing a physical reference system must necessarily lead 
to corrections compared to gauge fixing, because the physical reference system will communicate via grav- 
itational interaction with the original system under consideration. In the sense described above, the dust 
reference system creates only a minimal modification — it is the minimal extension of the original gravity 
- matter system that extracts the true degrees of freedom and allows for their physical evolution. 

The other difference is that instead of having constraints imposed on the phase space variables, G = 
Gj = 0, now the dynamics of the true degrees of freedom is subject to conservation laws H = Cj = 0. Thus, 
in solving (j4.5p we may prescribe arbitrary functions e{a), ej(fT) which play the role of the (constant in 
r-time) energy and momentum density, respectively. The substitutioiJ^ H = e, Gj = —ej will be crucial in 
what follows. In fact, in order to derive the second order equations of motion, H, Qjk in ()4.5p has to be solved 
for Il,P^^. Without the conservation laws, this would be impossible, since Il,P^^ enter the expressions for 
H, C, Gj in a non trivial way, i.e. solving for them would lead to algebraic equations of higher than fourth 
order. The substitution will also be crucial for the derivation of the effective Lagrangian, by the inverse 
Legendre transform, corresponding to H, see appendix [Cl 



'^^ For the explicit calculation of the Lie derivative it is important to note that 11 , P-''' are tensor densities of weight one in 
dust space. 

Note that we have included a cosmological constant term +2-^det(q))A in Cgco 

The letter e is chosen to indicate that these values are small, appropriate for test clocks and rods. In this way it can be 
guaranteed that the dust, although gravitationally coupled with the original system, will not alter the dynamics of the original 
system in an uncontrolled fashion. 



32 



4.2 Second Order (Lagrangian) Formulation 



In this section we will use the first order (Hamiltonian) equations of motion and derive the corresponding 
second order (Lagrangian) equations of motion for the configuration variables H and Qjk, respectively. We 
will sketch the main steps of these calculations in section [5]3j The reader who is just interested in the results 
should skip this section and go directly to section S3] where the final equations are summarised. 



4.3 Derivation of the Second Order Equations of Motion 



In this section we want to derive the second order equations of motion for H and Qjk, respectively. These 
second order equations will be functions of the configuration variables H, Qj^ and their corresponding veloc- 
ities 3,Qjk, respectively. This can be achieved by solving for the conjugate momenta 11 and P^^ in terms 
of their corresponding velocities H and Qjk via the equation of motion. The relation between the conjugate 
momenta and their velocities is given through the first order Hamiltonian equations which were displayed 
in the last section in equation (14. Sp . 



We begin with the matter equation for H. First, we have to take the time derivative of the first order 
equation for H given in equation (j4.5p . This yields 
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(4.8) 



As discussed in section \3.3\ the shift vector Nj := —Cj/H is a constant of motion since Cj = H = 0. 
Therefore for the Lie derivative with respect to the only non vanishing contribution is the one including 



(4.9) 



We will use this result later on, but for now we will work with the compact form of the Lie derivatives as 
written in equation (14. 8p . Solving for 11 in terms of H we get from equation ()4.5p 
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and thus have expressed 11 as a function of the velocity H. In order to stress that 11 has to be understood 
as a function of H, we have explicitly written the function's arguments in this section. Notice that strictly 
speaking 11 also appears in N = y'l + Q^'^CjCk/H^ and = -Q^kCk/H. However, Cj and H are treated 
as constants of motion as discussed before. The same applies to P^^ below. Next we insert this result into 
equation ()4.5p . obtaining 
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The final second order equation of motion for E can be derived by inserting equation (|4.1Up and (j4.1ip into 
equation ()4.8p . The result is 
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(4.12) 
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The same procedure has to be repeated for the gravitational equations now. Applying another time derivative 
to the first order equation of Qjk in equation (j4.5p yields 
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(4.13) 

In order to solve for P^^ in terms of Qjk we use the invcrsG of tliG tensor Gjkmn 

denoted by [G-^]jkmn 

defined in equation (j4.7|) . This results in 
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(4.14) 



Since the equation for P^^ in (j4.5p contains C which includes the geometry as well as the matter part of 
the Hamiltonian constraint (see equation (j3.56p for its explicit definition), it is a function of the variables 
Qjk,P^'',^ and n. This was different for 11 where its time derivative involved the matter momentum only. 
Thus, in order to express P^^ as a function of configuration variables and velocities, we use equation (j4.10p 
and (I4.14[) and replace the momenta occurring in P^'^. Rewriting Cgeo by means of the DeWitt bimetric 
Gjkmn we get 
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Using the relation in equation (j4.14p and the fact that G 

C,eo{Qjk,Qjk) = ^[^^^p[G-iF^™"(Q,fc- (£^Q).J(Q„„- (£^Q)„J +7d^ (4.16) 
For the matter part of the Hamiltonian constraint we obtain by means of equation (I4.10p 
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There are two other terms in P^^ which include the conjugate momenta P^^ . One is the first term on 
the right hand side of equation (|4.5p being quadratic in P^^ and the second is the Lie derivative of P^^ . 
Reinserting into those terms the relation shown in equation (j4.14p . we end up with the following expression 
for P^^ as a function of configuration and velocity variables: 
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Next we insert the expressions for P^'' and P^^ in equations (j4.14p and (j4.18p . respectively, into equation 
(j4.13p for Qjk- Keeping in mind that GjkmnQ"^"' = —^f^Qjk, we end up with 



Q 



jk 



1 1 mtrs F/^— 1 1 nuvw 



•la 



-limnvw i/^—l]turs 



[G- 



Qjk 



2VditQ ^ 2A ^ 



A 



2R 



jk 



NHk 



+ 2N{D,DkN) - -^^^ajkmuN"" 



Here we used that GjkmnlG^^]"^"-'''' = -Gjkmn[G~^y^'''' which foUows from {Gjkmn[G^^ 
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A straightforward, but tedious calculation shows that the second term on the right-hand side of the equation 
for (jjfc simplifies greatly: 
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Consequently, the final form for the second order equation of motion for Qjfc is given by 
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This finishes our derivation of the (general) second order equation of motion for E and Qjk, respectively. 



4.4 Summary of Second Order Equations of Motion 

The second order equations of motion for the (manifestly) gauge invariant matter scalar field E and the 
(manifestly) gauge invariant three metric Qjk have the following form: 
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and 
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The term C = Cgeo + Cmatter occurring in the equation for Qj^ has to be understood as a function of 
configuration and velocity variables and is explicitly given by 
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Apart from the fact that we have a dynamical lapse function given hj N = C/H, as well as a dynamical 
shift vector defined as Nj = —Cj/H, the only deviation from the standard Einstein equations in canonical 
form is the last term on the right-hand side in equation (j4.24p . This term, being quadratic in Nj and 
therefore quadratic in C,- = Cf^° + (^jnatter vanishes for instance when FRW spacetimes are considered. In 
our companion paper [15], we specialize these equations precisely to this context and show that the equations 
above reproduce the correct FRW equations. 



5 Physical Hamiltonian, Boundary Term and ADM Hamiltonian 

As long as the dust space S (and therefore X) has no boundary, H is functionally differentiable, which we 
always assumed so far. However, for more general topologies we are forced to consider boundary conditions. 
In this section we show how to deal with asymptotically flat boundary conditions for illustrative purposes. 
More general situations can be considered analogously. 

Recall |37] that asymptotically flat initial data are subject to the following boundary conditions 

Qabix) = 5., + + 0(r-2) , 

r 

r 

e(x) = ^ + 0{r-'), 

vr(x) = ^ + 0(r-3). (5.1) 

Here, x*^ is an asymptotic coordinate system, = 5abX°'x^ and Q denotes the angular dependence of the 
unit vector x"'/r (on the asymptotic sphere). The functions fab, /, g"''', g are assumed to be smooth. 
Moreover, fab is an even function under reflection at spatial infinity on the sphere, while g""^ is odd {f,g 
do not underly parity restrictions). Notice that these boundary conditions directly translate into analogous 
ones for the substitutions Qjk ^ Qab, P'^^ ^ P""^, ^ ^ H <-> vr, because switching from to 5 is just a 
diff eomor phism . 

The part of the differential of H that gives rise to a boundary term is 

'^boundary 

H = y d^a [iV^boundaryC + iV'' ^boundaryCj] , (5.2) 
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which coincides precisely with the boundary terms produced by the canonical Hamiltonia: 



-f^canonical 



Jx 



nc + n Cn 



(5.3) 



without dust. Here, the lapse and shift functions n, are Lagrange multiplier rather than dynamical quan- 
tities like . Therefore, the usual boundary terms [37] can be adopted, once the asymptotic behaviour 
of the dynamical lapse and shift functions N, have been determined, which, in turn, is completely spec- 
ified by Nj because = Q^^N^, N = y'l + Q)^NjNk. The scalar field contribution to C,Cj decays as 
thus it is sufficient to consider the geometrical contribution. Cj^° = —ID^Pj decays as and is 
even asymptotically. The term quadratic in P^^ and the term quadratic in the Christoffel symbols in C^'^° 
decays as while the term linear in the Cristoffel symbols decays as and is even. We conclude that 
(with H = ^JC^ — Q^'^CjCk) is asymptotically constant and even. The same is true for N, which 



H 



is anyway bounded from below by unity. The usual computation [22l [37] yields 
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+ I 7dit(Q)Q^'''^ 
2 / dS, NHpI , 



(5.4) 



where dSj = /rdQ , dQ the volume element of S'^. The prime is to indicate that, contrary to what 
the notation suggests, the terms shown are, a priori, not total differentials. In the usual formalism they 
are, because lapse and shift functions are Lagrange multipliers and do not depend on phase space. Here, 
however, they are dynamical and we must worry about the variation SNj. 

It turns out that the boundary conditions need to be refined in order to make H differentiable. We 
will not analyze the most general boundary conditions in this paper, but just make a specific choice that 
will be sufficient for our purposes. Namely, we will impose in addition that Cj^° = —2DkP^ decays as 



3+e 



1/r 

also H = 
decays as 



(e > 0) at spatial infinity. Then also Cj falls off as l/r'^^'^. Since C decays like 1/r'^, it follows that 
— Q^^CjCk decays as whence Nj = —Cj/H decays as l/r*^. This implies that 5Nj 



/r^ Thus, 5N = [Qi^NkSNj + NjNkSQ^^ /2]/N decays aslfH l/i 



2e 



It follows that 5B'{N) = and 



ti5B'{N) = kSEadm, Eadm = [ Vdet(Q)Q^'^ [-dSj T^k + Sk5r% 
reduces to the variation of the ADM energy. The correct Hamiltonian is given by 

H = Eadm + ^ dV H{a) . 



(5.5) 



(5.6) 



It is reassuring that in the asymptotically fiat context, we have automatically a boundary term in the Hamil- 
tonian, which is just the ADM energy. The additional bulk term comes from the dust energy density and 
does not vanish on the constraint surface. 



2^With the substitutions S ^ X, Qjk ^ Qab, P^'° ^ p"', H ^ ^, 11 ^ tt implied. 
■^"We could more generally have imposed that 5Nj falls off like l/r*^ at spatial infinity. 

■^^The choice e = 1/2 seems to be appropriate in order to reproduce the asymptotic Schwarzschild decay for dynamical lapse 
and shift. However, this is not the case here, because we are automatically in a freely falling (dust) frame and r is eigentime. 
Hence, Qtt ~ —N'^ + Q-'''NjNk — 1, whatever choice for Nj is made, it is independent of r. 
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Before we close this section, let us also mention the physical Hamiltonian system under consideration 
there exists a Lagrangian from which derives by Legendre transformation on the phase space defined by the 
physical observables. Curiously, the corresponding effective action turns out to be non-local in dust space 
but local in dust time. It can be computed via a fixed point equation to any order in the spatial derivatives 
of the fields. Details can be found in appendix O 



6 Linear, Manifestly Gauge Invariant Perturbation Theory 

In section [J] we derived the (manifestly) gauge invariant second order equation of motion for the scalar field 
S and the three metric Qjk- Now we want to consider small perturbations around a given background whose 
corresponding quantities will be denoted by a bar, namely H, Q^^. The linear perturbations are then defined 
as SE := E — E and 5Qjk ■= Qjk — Qjki respectively. Note that these perturbations are also (manifestly) 
gauge invariant because they are defined as a difference of two gauge invariant quantities. Consequently, 
any power of 5E and 6Qjk will also be a (manifestly) gauge invariant quantity such that in the framework 
introduced in this paper gauge invariant perturbation theory up to arbitrary order is not only possible, but 
also straightforward. This is a definite strength of our approach compared to the traditional one, see section 
[7] for a detailed discussion. 

However, in this section we will restrict ourselves to linear (manifestly) gauge invariant perturbation theory. 
That is any function F{Qjk,E) will be expanded up to linear order in the perturbations 6E and 6Qjk- We 
denote by 5F the linear order term in the Taylor expansion of the expression F{Qjk,E) — F[Qjj.,E). A 
calculation of higher order terms will be the content of a future publication. Usually, in cosmological per- 
turbation theory one chooses an FRW background and considers small perturbations around it. Here we 
will derive the equations of motion for the linear perturbations assuming an arbitrary background. In our 
companion paper, we will specialise the results derived here to the case of an FRW background and show 
that we can reproduce the standard results as presented, e.g., in 0. The reader who is only interested in 
the final form of the perturbed equation of motions should go directly to section 16.31 where a summary of 
the results is provided. 



6.1 Derivation of the Equation of Motion for 

Let us go back to the second order equation of motion for H shown in equation (j4.23p . Since its perturbation 
involves several terms we will, for reasons of book keeping, split this equation into several parts and consider 
the perturbation of those parts separately. Comparing the equation of H with the one for Qij in equation 
(j4.24p . we realise that in both equations the first term on the right-hand side includes an identical term, 
namely the expression in the square brackets. Therefore it is convenient to derive its perturbation in a 
closed form such that the result can then also be used for the derivation of the equation of motion of 5Qij . 
Throughout this section we will repeatedly need the perturbation of \/det Q and its inverse, given by 



^VdetQ = ]^^<leiQQ"HQij 

5(^=i=) = -\^=^-^'5Q,j. (6.1) 
Vdetg' 2 ^det Q 

Next, considering the definition of the lapse functions given by = C/H and the definition H = 
y/C^ - Q'-^CiCj we obtain N = ^Jl + Q^iNiNj. Thus, as mentioned before, N is not an independent 
variable, but can be expressed in terms of Qij and the shift vector, which itself is a function of the elementary 
variables H, Qij. However, as shown in our companion paper [15], the perturbation of Nj is again a constant 
of motion, that is 5Nj = 0. For this reason it is convenient to work with 6Qjk, SE and SNj, although, strictly 
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speaking, 6Nj is not an independent perturbation. Keeping this in mind the perturbation of yields 
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(6.2) 



From the exphcit form of 5N and 6^/6etQ we can derive the following expressions which we will need below: 
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Here the derivative with respect to r is understood to act on everything to its right, including the perturba- 
tions. We also used that 5Nj is a constant of motion, so the term proportional to 5N.j does not contribute. 
In order to calculate the perturbation of the Lie derivative term occurring in the square brackets in eqn. 
(I4.23p . we compute 
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The Lie derivative term yields then 
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Similar to the r-derivative, the Lie derivative acts on all terms to its right, including the linear pertur- 
bations. For the moment we keep the Lie derivative with respect to 5N in the compact form above. At a 
later stage we will write down these terms explicitly and express them in terms of 5Qjk and 5Nj. Having 
calculated the variations of all components of the first square bracket in eqn. (j4.23p . we can now give the 
final result: 
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The last term in this equation can be written explicitly in terms of SQmn and 6Nm as 
N VdetQ \ r N 
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As the terms in equation (|6.6p are multiplied by (H— in equation (j4.23p . we also need the perturbation 
of the latter term. It is given by 
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dr N 



(6.8) 



Next we determine the perturbation of Q^^E f^ which enters the second term on the right-hand side of 
equation (|4.23p : 
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We will split the third term occurring on the right-hand side of equation (j4.23p into N'^ and the remaining 
part given by (AS + [Q^^]jE^k — Their respective perturbations are 
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Finally, we have to calculate the perturbation of the last three terms in equation (j4.23p . involving Lie 
derivatives of E and S, respectively. We obtain 
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Adding up all the contributions and factoring out the linear perturbations 6Qjk and SNj, we can rewrite 
the second order EOM as 
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Here we used that the last term occurring in equation (j4.23p . which involves the Lie derivative with respect 
to (5A, can again be expressed in terms of the perturbations 5Qjk and 5Nj, explicitly given by 
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Analogously, the last term on the right hand side of equation (j6.6p . which also involves a Lie derivative with 
respect to SN, can be expressed as 
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Moreover occurring in equation (j6.8|) is given by 

5Nj. (6.16) 

We would like to emphasise again that the partial and Lie derivatives act on everything to their right and 
therefore also on the perturbations. That is the reason why for instance (5H and 5Qjk have not been factored 
out separately. This finishes our derivation of the second order equation of motion for the scalar field 
perturbation JS. In the next section we will discuss the corresponding equation for the metric perturbation 
6Qjk- 



Q 



5Qjk + 



6.2 Derivation of the Equation of Motion for 6Qjk 

Similar to the derivation of the second order equation of motion for 6'E we will also split the equation for Qjfc 
in equation (j4.24p into several terms whose perturbations are then considered separately. More precisely, 
we will split the equation into seven separate terms, with the three terms involving Lie derivatives in the 
last line counted as one. Starting with the first term on the right-hand side, we recall that the perturbation 
of the sum of terms in the square brackets has already been computed during the derivation of the equation 
for (5H. Thus, we can take over those results, as shown in equation (j6.6p . The perturbation of the term 
{Qjk - {C^Q)jk) is given by 
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Consequently, we have all ingredients needed for the perturbation of the first term. However, we will not 
present the final expression in the main text since it is quite lengthly. Nevertheless, the interested reader 
can find the explicit form in appendix [G] in equation (j6.17|) . The perturbation of the second term on the 
right-hand side yields 
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The perturbation of the third term occurring on the right-hand side requires a bit more work, because it 
involves C = C^^° + (7™'^**'^''. Thus we will perform this calculation in two steps. First we will ignore the 
explicit form of 5C in terms of SQjk, (5^ and SNj. The resulting expression looks like 
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Due to its length, the exphcit calculation of 5C can be found in appendix El However, when actually 
inserting the expression of 5C into equation (j6.20p . some of the terms in the expression for 5C cancel with 
existing terms in equation (j6.20p . As a result, the final expression for the perturbation of the third term in 
equation (j6.2ip is slightly less involved. It is given by 
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Proceeding with the next term on the right-hand side of equation (j4.24p . we obtain for its perturbation the 
following result: 
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For this calculation we used the fact that the perturbation of the Ricci tensor can be expressed in terms of 
the perturbations SQjk- The explicit relation reads 
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The next term in equation (j4.24p includes covariant derivatives. Therefore we will have to consider the 
perturbation of the Christoffel symbols T^. These can again be written in terms of metric perturbations as 



shown below 

Using this we end up with 

5{2NDjDkN) 



STTk 



^ ^{DjSQnk) + {Dk5Qjn) - (DnSQjk 



AN(D^DkN) - 2N D^Dk 



1 N N 



L2 ]y' 



+ 
+ 



m(DjDkN) ]+2N^DjDk 



INN 
2^ 

SNn, 



6Q, 



(6.24) 



(6.25) 



w 

2N{DnN)D^kT''] SQj),n. + \-N^ {DnN)DmQ"''] 6Qjk. 



43 



Note that covariant derivatives surrounded by round bracket such as {DjDj....) act on the elements inside the 
round brackets only. By contrast, covariant derivatives not surrounded by round brackets act on everything 
to their right, including also the perturbations. 

Next we deal with the three terms on the right-hand side of equation (j4.24p that include Lie derivatives 
with respect to 5N . Those Lie derivatives are again functions of 5Qjk, SNj, and partial derivatives thereof, 
because we have SN'^ = —Q Q N s^Qm + Q 5Nn- In this section, however, we will work with the 
following compact form only: 
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In the next section, we will rewrite the second order equation of motion in a more concise form, using 
coefficient functions. That will allow us to include the Lie derivative term with respect to 5N in the 
following, more explicit form: 
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Another term which includes Lie derivatives with respect to 5N appeared previously as part of equation 
(j6.19p . Performing the Lie derivative also in this case, we end up with 
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The other two terms that include Lie derivatives with respect to 5N are parts of the perturbation of the 
first term. Their explicit expressions are given in equation (j6.7p and equation ()6.18p . respectively. Note that 
these terms are written out explicitly in appendix [Gl where the final form of the perturbation of the first 
term is calculated. 

Finally, we consider the last remaining term from equation (|4.24p . It involves the Hamiltonian density H, 
which we found out to be a constant of motion in section (3.4). In our companion paper [15] we show that 
also 5Nj and 6H are constants of motion. Therefore, in complete analogy to the case of 5Nj, we will factor 
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out 5H. We thus obtain 
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Now we have finally derived all the individual parts that are needed in order to write down the equation of 
motion for SQjk- However, by just looking at the various individual terms, it is clear that they are already 
considerably more complicated than for the corresponding case of the matter equation of motion for 6'B. 
Nevertheless, we decided to present the final equation in detailed form on the next page, in particular to 
convey a sense of how much more involved the geometrical part of the perturbed equations is compared to 
the matter part. In the next section we will rewrite both equations of motion, the one for (5H as well as the 
one for 5Qjk, in a more transparent form where all the complicated background coefficients in front of the 
perturbations are hidden in certain coefficient functions. In our companion paper, we will then specialise 
those coefficient functions to the case of an FRW spacetime and show that the general equations derived 
in the last two subsection are (up to small correction caused by our dust clock) in agreement with the 
well-known perturbation equations as discussed, e.g., in [7]. Note that we still kept the compact form of 
the Lie derivative with respect to SN in equation (j6.3ip . because we wanted to present this equation on one 
page only. 
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6.3 Summary of the Equations of Motion for 5S and 5Qjk 

In the last two sections we derived the second order equations of motion for 5E and SQjk- The results of 
our calculations can be found in equation (j6.13p and equation (j6.3ip . respectively. However, these equations 
are quite complex and not very transparent in that form. For this reason, we want to rewrite them in 
a form where we can still recognize their general form but where they take a much simpler form. We 
will hide the precise details of the various background quantities that occur as coefficients in front of the 
linear perturbations in certain coefficient functions that will be introduced below. These coefficients will 
be operator valued since they also involve objects such as partial or Lie derivatives, as can be seen from 
eqns. ()6.13p and ()6.3ip . As explained before, apart from the elementary perturbations dQjk and 6E, the 
second order equations of motion contain also the perturbation of the shift vector 5Nj and the (physical) 
Hamiltonian density 5H. Both are functions of SQjk and 53, and therefore not independent perturbations. 
However, it turns out that these two quantities are constants of motion, so it is convenient to keep them in 
the equations. 

Starting with the second order equation for 5H, its general structure is of the kind 



(6.32) 



where the coefficients are given by 



9r2 



N^dT 



N {^/detQ) ^ N 



N y/detQ VdetQ^ ^ 



C-. 



VdetQ 



-Q 



N 



VdetQ 



N J det Q),j 



d 



-r-f2 / . r-^mn, d 



(6.33) 



and 



a 



i- 



d_ ^ V det Q -Qjmj^'- 



-p^jm-p^kn-r-: TTz^ / - — \ -I -p^jm-p^kn— r N 



+ 



\N (VdetQ)' 



+ 



N 



C- 



VdetQ 



N VdetQ VdetQ ^ ^ N 

+ (Q =,n) 

N 



N 1 f-^jk N^n'' 



VdetQ 



N 



8 /l^Q.^_^)^^detQ 



N 



(6.34) 



(A VdetQ), 



Note that all partial and Lie derivatives act on all terms to their right, including the linear perturbations. 
This is also the reason why terms as for instance H or Qjk do not occur in the simple form of equation 
(1632)1 . 
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The third coefficient is given by 
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As was to be expected, the corresponding equation for the perturbation of the three metric SQjk includes 
more than just three terms. It is of the form 
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The various coefficients introduced in the equation above are given as follows: 
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Here we use the notation that covariant derivatives surrounded by round bracket such as {DjD]....) act on 
the elements inside the round brackets only. By contrast, covariant derivatives not surrounded by round 
brackets act on everything to their right, including also the perturbations. The coefficient for 5E can be 
explicitly written as 



(6.38) 

For the coefficient belonging to the linear perturbation of the Hamiltonian density 5H we get: 
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The last two coefficients are the ones for SQmn and SNm- These are the most complicated ones for the 
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second order equation of motion of 6Qjk- We will list them below: 
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Note that the coefficient for 6Qmn in equation (j6.3ip and the one above are different due to the presence 
of the last two lines in the equation above. The reason for this is that now we used the explicit expression 
for the Lie derivatives with respect to 6N, which were derived in equation (j6.28p and (j6.27p and lead to 
additional terms in 6Qmn and 5Nm- 

Finally, we present the coefficient for 5Nfn. Similarly to the case of SQnm, we also get additional terms 
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coming from the Lie derivatives in the last hne of equation (|6.31|) . 
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Although the form of the perturbed equations is quite complicated, they simplify drastically for special 
backgrounds of interest. For the case of FRW, for instance, all terms proportional to Nj vanish, since 
N j = —Cj/H = for FRW. This is due to the geometry and matter parts of the diffeomorphism constraint 
vanishing both separately in that case. Furthermore, all terms in the coefficients that contain spatial deriva- 
tives applied to background quantities vanish also. Other backgrounds where considerable simplification 
will occur include Schwarzschild spacetime. 



7 Comparison with Other Frameworks 

We now proceed to compare our work with other approaches to general relativistic perturbation theory found 
in the literature. In the following we will restrict ourselves to discussing works that treat perturbation theory 
around general backgrounds. Approaches which deal exclusively with cosmological perturbation theory will 
be looked at in our second paper, specifically dedicated to that topic. 

The central point of comparison is how gauge-invariance is handled in the various approaches. As that 
notion often acquires different meanings, especially in the context of general-relativistic perturbation theory, 
it seems prudent to recall the precise mathematical setting underlying most works. It was developed by 
Sachs [27] and Stewart and Walker [28j, and recently given a very general and elegant formulation by Bruni, 
Sonego and collaborators [29l|30ll3l]. The starting point consists of two spacetime manifolds Mq and M, 
which represent the background spacetime around which one perturbs and the actual physical spacetime, 
respectively. It is important to keep in mind that Mq is only an artificial construct. Perturbations of 
geometric quantities are then defined by comparing their values in Mq and M, respectively. This procedure 
is highly ambiguous, however, in that there is a great freedom inherent in the choice of points where one 
compares background and "real" quantities to define the perturbations. Note that this freedom is in addition 
to the usual coordinate gauge freedom in general relativity. For that reason, Sachs termed it gauge invariance 
of the second type |27j . Making such a choice, which mathematically amounts to choosing a so-called point 
identification map between Mq and M, is therefore nothing but a choice of gauge. Correspondingly, gauge- 
invariant perturbations are those quantities whose values do not depend on the choice of point identification 
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map. The actual condition for a perturbed quantity to be gauge-invariant to first order in this sense was 
aheady derived in ^27j, fully proved in [28j and finally generalized to arbitrary order n in |29j . The result 
can be succinctly summarised as follows: a geometric quantity T - such as a tensor - is invariant to order n 
iff all its perturbations to order n — 1 are either vanishing, (spacetime) constant scalars or a combination of 
Kronecker deltas. This result is often known as the Stewart&; Walker lemma. Clearly, the only case of actual 
interest is the first. An example is given by curvature tensors in linear perturbation theory around Minkowski 
space. As they vanish in the background, they are gauge-invariant to first order. These insights have been 
the backbone of most attempts to construct gauge-invariant quantities to various orders in perturbation 
theory or even non-perturbatively. We will now discuss two of them. 

In a series of papers [21 El Sj Nakamura has used these principles to develop formulas for gauge-invariant 
quantities to second and third order around an arbitrary background. They encompass the invariant parts of 
various metric and curvature, as well as matter perturbations. These general formulas are, however, implicit 
only to the extent that Nakamura derives them from the assumption that the corresponding linear order 
perturbations can be decomposed into gauge-invariant and gauge-variant parts. Consequently, while the 
construction is, in principle, valid for arbitrary backgrounds, in practice only backgrounds with sufficient 
symmetries to perform that split at linear order explicitly can be used. Luckily, that applies, of course, to 
several cases of great interest, such as cosmological and spherically symmetric backgrounds. The latter case 
is explicitly treated in [4]. 

A distinctly different approach - let us call it the EB approach - is based on seminal work by Ellis and 
Bruni [32] , which has since inspired a multitude of other works \33\ [Ml [35] . Although the discussions in 
these papers is geared towards applications in cosmology, the framework itself can be applied to arbitrary 
spacetimes, in principle, which is why we decided to discuss it here. The basic idea is to use a 1 + 3- 
approach by employing covariant quantities connected to a family of flow lines or " fundamental observers" . 
The prime reason is that these quantities are much more closely related to what one actually measures in 
astrophysics. Furthermore, by a simple application of the Stewart- Walker lemma, they are automatically 
gauge-invariant if the corresponding quantities in the background spacetime vanish. Unlike in the more 
common metric-based formalism, these quantities are defined in the physical, perturbed spacetime. As a 
result, they are fully non-perturbative and in that sense their gauge-invariance extends to all orders. The 
connection to the standard perturbative approach based on perturbations in the background spacetime can 
be made by suitably expanding the physical quantities to the desired order, see |33J- This approach enjoys 
a clear geometric and physical interpretation of the quantities used, as well as the advantage of basing 
perturbation theory on non-perturbative variables. 

Comparing the works mentioned so far to ours, a first obvious difference is that we work in a canonical 
setting versus the covariant setting used by the others. The motivation is first that gauge issues become 
particularly clear in the canonical picture and second our view towards quantization. The more important 
difference, however, is our use of dust as a dynamically implemented coordinate system. Our dust clocks 
serve a twofold purpose: on the one hand they enable us to construct background observables and therefore 
to solve the standard gauge problem in general relativity. On the other hand they also serve as the point 
identification map and thus eliminate the gauge freedom of "second type". In that sense they represent a 
logical extension to perturbation theory of the initial conceptual idea by Brown and Kuchar [12] to use dust 
as a physical and therefore preferred coordinate system. 

We should also point out that while our framework employs the metric and its perturbations as funda- 
mental variables, we could equally well use the dust clocks to build a gauge-invariant perturbation theory 
based on the same variables used in the EB approach. In fact, it seems worthwhile to look a bit more closely 
at the relationship between the two approaches. Both are non-perturbative in the following sense: they con- 
struct quantities which are gauge-invariant (with respect to gauge transformations of the second type). Only 
then perturbation theory is applied, which means that gauge-invariance is then automatically guaranteed 
in each order of perturbation. The difference arises when one looks at the role of gauge transformations of 
the first kind. The EB approach uses idealized observers that are comoving with the physical matter in the 
model. Thus the theory is not deparametrized and gauge freedom with respect to the background spacetime 
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remains, as illustrated by the presence of constraints as part of the equations of motion. In our case, the 
observers are represented by the dust, a component added to the physical matter content of the theory. 
They are thus dynamically included in the theory via the dust contribution to the Lagrangian, in addition 
to all the other matter. This allows for a complete deparametrization of the non-dust system with respect 
to the dust. Time evolution for this subsystem becomes unconstrained and a physical Hamiltonian emerges. 
The price to pay for this is that the dust contributes to the energy-momentum tensor of the deparametrized 
system, the size of which is small, however. One might well argue that, for practical purposes at the classical 
level, the choice between the two approaches is a matter of taste. Our approach, however, offers a clear 
advantage if one is interested in quantization, eventually. All programs aiming at a quantization of gravity 
that have been pursued, so far, use the metric as a fundamental building block. It is not obvious to us how 
to attempt a quantization based on the covariant variables used in the EB approach. 

Another argument in favour of our framework is the following: A test observer which by definition does 
not have any impact is only a mathematical idealisation. Physically much more realistic is a dynamically 
coupled observer fluid like the dust considered here which in particular takes into account the gravitational 
backreaction. 

To summarise: the crucial difference between our approach and the others discussed here is that the latter 
deal only with gauge freedom of the second type. This can be seen from the fact that they use background 
variables which are not gauge-invariant and is evidenced, e.g., by the presence of constraint equations in 
addition to evolution equations. Our treatment, by contrast, deals with all variables at all orders in a unified 
manner. Furthermore, in our opinion the framework developed here allows for a much more straightforward 
implementation at higher orders. Another advantage is that it allows to treat arbitrary backgrounds in 
practice, without the high degree of symmetry necessary for approaches based on the Stewart- Walker lemma 
to work. Recall that the latter require finding non-trivial quantities that vanish in the background manifold. 
Only in symmetric backgrounds such as homogeneous spacetimes is that a fairly tractable problem. 

Finally, we should briefly discuss the recent work in [25], which is close to ours in terms of motivation 
and conceptual underpinnings. The authors there also use the general gauge invariant framework of 
however, with two differences: First of all, they do not use dust matter to achieve gauge invariant completions 
of geometry and matter variables. This prevents them from bringing the constraints into a deparametrised 
form [9] and thus there is no time independent physical Hamiltonian. Secondly, while they can develop 
higher order cosmological perturbation theory, their perturbations of gauge invariant quantities are still 
expanded in terms of the perturbations of the the gauge variant quantities which is what we never do. 
Therefore the basic perturbation variables are different in the two schemes: In our scheme we never care 
how our gauge invariant variables are assembled from gauge variant ones, they and their perturbations are 
fundamental for us and n'th order quantities are n'th order expressions in those. In contrast, in [25j n'th 
order means n'th order in the gauge variant quantities. In particular, the n'th order perturbed variables 
are only invariant with respect to the n'th order constraints up to terms of higher order. In contrast, our 
perturbed variables are always first order and always fully gauge invariant, it is only in the Hamiltonian that 
higher orders of gauge invariant variables appear. It would no doubt be fruitful to translate the schemes 
into each other and to see which differences and similarities arise. 

8 Conclusions and Open Questions 

This is a long and technically involved paper. The reader rightfully will ask why one should dive into its 
details and what exactly is novel as compared to the existing literature. The following remarks are in order: 

1. Non perturbative gauge invariance 

To the best of our knowledge, there exists no generally accepted notion of gauge invariance at n-th 
order of perturbation theory in General Relativity. Moreover, at each order of perturbation theory 
one has to repeat the analysis for how to preserve gauge invariance to the given order. Given those 
difficulties, it is natural to try to invent a scheme which separates the issue of gauge invariance from 
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the perturbation theory. Hence, one must treat gauge invariance non perturbatively. This is exactly 
what we managed to do in this paper. 

Thus, one only deals with the exact observables of the theory. All the equations of the theory have to 
be written in terms of those gauge invariant quantities. Given such a gauge invariant function O on 
the full phase space, we evaluate it on a certain background (data) which is an exact solution to our 
equations of motion and get a certain value Oq. The perturbation of O is then defined as 50 = O — Oq. 
We never care to expand 50 in terms of the perturbations of the gauge variant degrees of freedom 
(although we could). However, all the equations are expanded directly in terms of the perturbations 
of those physical observables. 

2. Material reference system 

In General Relativity it is well known that in order to meaningfully talk about the Einstein equations 
and to have them describe something observable or measurable, one has to suppose that spacetime is 
inhabited by (geodesic) test observers. By definition, a test observer has no effect whatsoever on the 
system. This is of course mathematically convenient but physically worrysome because a test observer 
is a mathematical idealisation. Every real observer interacts at least gravitationally and does leave 
its fingerprint on the system. One of the achievements of the seminal work [12] of Brown & Kuchaf , 
which in our mind is insufficiently appreciated in the literature, is to have overcome this shortcoming. 
The authors of [12] have identified a generally covariant Lagrangian which comes as close as possible 
to describing a non - self interacting, perfect and geodesically moving fluid that fills out spacetime 
(congruence). It does leave its fingerprint on the system and thus is physically more realistic than the 
test observer fluid. 

In this paper we have driven the work of [12j to its logical frontier and have asked the question 
whether the dust when added to the geometry - matter system really accomplishes the goal of keeping 
the approximate validity of the usual interpretation of the Einstein equations. We have verified that 
it does which in our mind is an intriguing result. 

3. Solving the Problem of Time 

Since General Relativity is a generally covariant or reparametrisation invariant theory, it is not 
equipped with a natural Hamiltonian. Rather, the "dynamics" of the non observables is described 
by a linear combination of constraints which really generate gauge transformations rather than phys- 
ical evolution. Observable quantities are gauge invariant and therefore do not evolve with respect to 
the "gauge dynamics". Therefore it is conceptually unclear what to do with those observables. The 
achievement of |10^ [TT] is to have invented a scheme that in principle unfreezes the observables from 
their non motion. However, in general that physical motion is far from uniquely selected, there are 
in general infinitely many such physical notions of time and none of them is preferred. Moreover, the 
associated Hamiltonians are generically neither preserved nor positive or at least bounded from below. 

When combining the frameworks of [12] and [10^ [TT] we find the remarkable result that there is 
a preferred Hamiltonian which is manifestly positive, not explicitly dependent on physical time and 
gauge invariant. It maps a conceptually complicated gauge systems into the safe realm of a conservative 
Hamiltonian system. The physical Hamiltonian drives the evolution of the physical observables. It 
reproduces the Einstein equations for the gauge invariant observables up to corrections which describe 
the influence of the dust. 

4. Counting of the Physical Degrees of Freedom 

The price to pay is that one has to assume the existence of the dust as an additional matter species 
next to those of the standard model^. It would maybe be more desirable to have matter species of 

Curiously, what we have done in this paper bears some resemblance with the Stueckelberg formalism. There one also adds 
additional matter to Maxwell theory. One can then make the longitudinal mode gauge invariant and thus finds a theory with 
one more degree of freedom. This is one way to arrive at the Proca theory and more generally at massive vector boson theories 
via the Higgs mechanism. 
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the standard model or geometry modes playing the role of a dynamical test observeijfl . In principle 
this is possible, however, the resulting formalism is much more complicated and it does not lead 
to deparametrisation. Thus a conserved physical Hamiltonian would then not be available and the 
equations of motion would become intractable. 

It is true that the dust variables disappear in the final description of the observables which are com- 
plicated aggregates built out of all fields. However, the theory has fundamentally four more physical 
degrees of freedom than without dust and that might eventually rule out our theory if those additional 
degrees of freedom are not observed. 

The truly remarkable feature of the dust is that it replaces the initial value constraints of General 
Relativity which are reponsible for having less physical degrees of freedom than one would naively 
expect, by four conservation laws. That is, in any given solution of our equations of motion, the 
physical observables must physically evolve with respect to each other in such a way that the conserved 
quantities do not change. This effectively acts like a constraint and therefore reduces the number of 
independently evolving observables by four, in agreement with the counting of the degrees of freedom 
without dust. Thus, at least as long as the value of those conserved quantities is sufficiently small, we 
will not be able to see those additional degrees of freedom. It is this fact which makes it possible that 
one effectively does not see more degrees of freedom than in the standard treatment. 

As a final objection against our formalism one might raise the fact that the dust contributes with the 
wrong sign to the matter energy momentum tensor. However, the formalism not only forces us to do 
this as we would otherwise have a negative definite Hamiltonian, moreover, as already remarked in 
the introduction, it is completely acceptable since here we talk about the energy momentum tensor of 
non observables. The energy momentum tensor of the observables in the final Einstein equations does 
satisfy the usual energy conditions. 

5. Complexity of the Equations of Motion 

Since General Relativity is a highly non linear, complicated self - interacting theory, experience from 
much less complicated integrable systems suggests that its Invariants, that is, the gauge invariant 
observables, satisfy a tremendously complicated Poisson algebra and that the equations of motion are 
intractable. Surprisingly, this is not at all the case. The observable algebra is almost as simple as the 
algebra of non observables and the equations of motion can be solved almost as easily as in the usual 
gauge variant formalism. Key to that is the presence of the already mentioned conserved current. 

In this first paper we have developed the general gauge invariant formalism and linear perturbation theory 
about general backgrounds. In the companion paper we apply the general results to fiat and FRW back- 
grounds and find agreement with usual linear perturbation theory for linearly invariant observables. This 
is a first consistency test that our theory has passed. Thus we hope to have convinced the reader that the 
present framework has conceptual advantages over previous ones and that it is nonetheless technically not 
much more complicated. Actually, the pay - off for having a manifestly gauge invariant approach will really 
come in at higher order where we believe that our equations of motion will be simpler. 

There are many lines of investigations that one can follow from here. An obvious one, the specialization 
to the all-important case of an FRW background is presented in a companion paper as already mentioned. 
Investigating perturbations around backgrounds of astrophysical interest, such as Schwarzschild spacetime, 
is also valuable. Again, for all these cases it should also prove very interesting to go beyond the linear to 
higher orders. Predictions from second-order perturbation theory, e.g. the issue of non-Gaussianity in cos- 
mological perturbations, are the topic of current research and could soon be testable by future experiments 
such as PLANCK, see [38]. On a more technical and conceptual level, our framework might prove useful 
to settle the issue of under what conditions general-relativistic perturbation theory is consistent and stable. 

'^'^This would be similar to technicolour theories which declare the Higgs scalar field not as an independent degree of freedom 
but as a compound object built from the bosons of the electroweak theory. Here one would build four independent scalars e.g. 
from the geometry field. 
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Finally, with a view towards facilitating the search for physically relevant predictions from approaches to 
quantum gravity, a quantization of our gauge-invariant formulation of general relativity, together with the 
development of the corresponding perturbation theory at the quantum level strikes us as a highly desirable 
goal. See [23] for first steps. 
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A Second Class Constraints of the Brown — Kuchar Theory 

In this section we provide the calculational details of the constraint analysis of the Brown - Kuchar theory 
discussed in section [2j In particular we want to show that no tertiary constraints arise. Our starting point 
is equation ()2.25p which we display once again below for the convenience of the reader 

Z,t — {-f^primary ) — C 



,tot 



Za,t — {-^primary ) Pa} — 
Z,t = {i^primary,/} = 5[ j=^^ + ^/d^){q-'UaUb + I)] =: C 

2 ^det(g)/)^ 

Zf, = {//primary, P} = -f^'P " npy^d^q^'UaS^^ + PS^^ 

Zj,t = {Hprirne.ry,Pj+WjP} = f,,P-{n''-'^^^y^p^q''''Ub)PWj,a (A.l) 



The last two equations involve the Lagrange multipliers and /Xj and can be solved for them. In contrast 
we observe that the first three equation are independent of the Lagrange multipliers, they are secondary 
constraints. We will now proceed with the constraint analysis and show that when the Poisson brackets 
between the primary Hamiltonian .ff primary and the secondary constraints are consideres no new constraints 
are generated. Recall that the primary Hamiltonian density was given by 

Vimary = fl^ Zj + flZ + flj + UZ + Za + u' c'°' + u'^c'"^ (A.2) 
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whereby the single constraints are shown below 

gtot _ ^eo _|_ ^matter _|_ 
gtot _ ^eo _|_ ^matter _|_ 



1. vr2 



c° = l[^^=- + Vd^p(g«''C/aC/6 + l)] 

Xc^-tter ^ ^^^^ 

c° = P[r„-w,.5y (A.3) 

We begin with the calculation of the Poisson bracket of -ff primary and the smeared constraint 

:= j (A.4) 



c 



{i?prima.y, c-*°*(n)} = j d^j d^y({^x\y)Z^{y), n\xy:\x)} + {^l{y)Z{y),n%x)c':\x)} 

+{l^^{y)Z\y),n\x)c'^\x)} + {u{y)z{y),n''{x)c'^\x)} + {i.^(y)2;;.(y), n«(x)4°*(x)} 
+{n'{y)c'°\y),n\x)c':\x)} + {n'^(y)4°*(y), (A.5) 

For the single Poisson brackets that occur in the equation above we obtain the following result: 

j d\ j d^y{ii^{y)Zj{y),n%x)c'':\x)} = j d\n''Pti^Wj,a (A.6) 

j d^x J d^y{fi{y)Z{y),n%xy:\x)} = 

J d^x J dMt^jiy)Z'iy),n''ix)c':\x)} = - j d^xPt,,S%n^ 

I d^x j d^y{u{y)z{y),n-{x)c':\x)} = 

d^x J dM^\y)z,{y),n-{x)ax)} = c^-\u) 

j d\ j d^y{n'{y)c'°\y),n%x)c'':\x)} = c'°\jCnn') + c{Cn' p) + J d^xpn' n''^/d^q''^Ub{Wj,aSi - Wj,cS^,) 

I' d'x J d'y{n'\y)cl-\y),n%xy:\x)} = ciC^n') - J d'x{n"'n' - n''n"')PS%Wj,, 



I 
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Consequently we can rewrite equation (|A.5p as 

{^primary, ^°*(n)} = c'°\C^n') + c(£^p) + j xpnn^ y^d^)q''^Uk{Wj^aSi - Wj^.Si) 

- J d?x{n'\'' - n"'n^)PWj^aSi^ 

+Wj^a ( - Pn''S{ + pn y^d^q'^UbS^, + f,^ P) ) (A.7) 



Hence, the result above involves the Langrange multipliers pj and p^ and can be solved for them such that 
no new constraints arise from c*°*. Proceeding with c*"* whereby the smeared constraint is given by 



ct°t(n) := J d^xn{x)c'"\x) (A.8) 

Thus we get 

{i^pnmary,c*°*(n)} = Jd^xJ d^y (^{fi^ (y) Z, (y) , n{x)c'°' (x)} + {p{y) Z (y) , n{x)c{x)} 

+{f,,{y)Z\y),n{x)c'°\x)} + My)z{y),n{x)c'°'{x)} + {u\y)zt{y),n{x)c'^^ 
+{n'{y)c'°'{y),n{x)c'°\x)} + {n"'(y)cr (y), n(x)c*°'(x)}) (A.9) 

For the single Poisson brackets that occur in the equation above we obtain the following result 



d^x j dMf^'iy)Zjiy),n{x)c'"\x)} = - j d^x^det{q)npfi^q^'UkWj,a (A.IO) 

d'x I dMKy)Ziy),n{x)c"'\x)} = I d^xp^~c 

d^x J d^y{pj{y)Z^{y),n{x)c'"\x)} = J d\^i,pn^/d^)q^''Ui,Si, 

d'x I d'yMy)ziy),n{xy°'{x)} = c'°\u) 

d^x I dS{iy'{y)zh{y),n{xy°\x)} = 

J d^x J d:Mn'{y)c'''\y),n{x)c'°'{x)} = ^°\q-^[n'dn - ndn']) 

d^x [ dh{n'\y)cf\y),n{x)c'°\x)] = -c*°*(£^,n) - c(£^,/>) 



d''xpnn"'^/d^q'"'Ub{Wj,aSi - H^i,c5f„) 
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Reinserting these results into equation (|A.9p yields 



j d^x^/d^)npq''^U,{Wj,,[n"'Sl - p^ - P^^d^q^^U^S^) 



S^,{pj - W.^au"' + P^q'^-UaW,, 



j d''x^jA^npq'''Ub{Wj,c{n"'S^^, - p^ - ^7dit(^g'^^C/rf5f,) 



S'A^'j - W,,an'^ + ^q^'UdW.^e)) (A.ll) 

Here we used in the last step that the last integral in the line before the last line one vanishes, because 
Wj^cS'',e — Wj^eS'lc is antisymmetric in e, c and multiplied by q^'^q'^^Ui^Ud which is symmetric in the indeces 
c, e. 

These are again the equation involving the Lagrange multipliers that we have seen before in the calculations 
for c*°*. 

Finally, let us consider the Poisson bracket of ^^primary and the secondary constraint c whose smeared 
version is given by 

c(n) := j d^xu{x)c{x) (A.12) 
where u is an appropiate smearing function. We obtain 

{-f^primary,c('u)} = j '^^^ j '^^V {{f^^ iv) ^ j iv) ^ + ^ iv) ^ 

+{f^j{y)Z^iy)^ u{x)c{x)} + {u{y)z{y), u{x)c{x)} + {i'^{y)zb{y),u{x)c{x)} 
+{n'{y)c'''\y),u{x)c{x)} + {n"{y)cl'''{y),u{x)c{x)}) (A.13) 

In this case we do not need to compute all the individual Poisson bracket in order to convince ourselves that 
no constraints arise, because the Poisson bracket of Z{p) and c{u) yields 

d^x [ d^y{p{y)Z{y),u{x)c(x)} = [ d^xpu "^Z^ (A.14) 
J J p'^y'det(g) 

which is a new term invloving the Lagrange multiplier p. Thus, we can solve the equation {-f^primaryj c{u)} = 
for p. 

It follows that no new terms are produced not involving p^ , pj, p in this second iteration step. Consequently, 
the full set of (primary and secondary) constraints is given by c*°* , c*°* , c, Zj ,Z^,Z, Za and z and it remains 
to classify them into first and second class. Obviously, 

{Z\x),Zk{y)] = P6i6{x,y), 

{Z{x),ciy)} = -^=5ix,y), (A.15) 
Vdet(g) 
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does not vanish on the constraint surface defined by the final set of constraints, hence they are second class 
constraints. Since n appears only linearly in c and does not appear at all, it follows that z, Za are first 
class. 

Let us consider the linear combination 

C*°* = I p,a + I' Wj^a + PTa + Pj + p n,, + Va + Ca 

= C*°* + Z p^a + Wj^a + Zj + Zn,a + Za (A. 16) 

where 

Ca = cr + C""'' (A.17) 

is the non-dust contribution to the spatial diffeomorphism constraint cjf*. Since all constraints are scalar or 
covector densities of weight one and cjf* is the generator of spatial diffeomorphisms, it follows that c*°* is 
first class. Finally, we consider as an Ansatz the linear combination 

gtot ^ ^tot ^ ^, ^ ^, + aZ , (A. 18) 

and determine the phase space functions such that c*°* has vanishing Poisson brackets with 

Zj, Z^ , Z up to terms proportional to Zj, Z^ , Z. 
We have 

{c*°*(x),Z,(y)} = {c*°*(x),Z,(y)}+afc(x){Z'=(x),Z,(y)} 

= {c'°\x),Z,{y)} + a,{x)P{x)5\x,y). (A.19) 

where we used equation ()A.15p in the last line. Solving this equation for aj we end up with 

a,{x) = - j d'y^{c'°\x),Z,{y)} = (^^^/d:^) pq'^'U.W,,^ (x) (A.20) 

which is a sensible expression since {c*°*(x), Zfc(y)} ~ 6^{x,y). For the Poisson bracket involving Z^ we get 

{c'"\x),Z^iy)} = {c'"\x),Z^{y)} + a''{x){Zkix),P{y)} 

= {c'"\x),Zj{y)}-a^{x)Pix)d^ix,y) (A.21) 

such that this equation can be solved for explicitly given by 

a^{x) = J d3y_L_{ctot(^)^ z^^y)} = (^^pV^}'l'"'UtS^^c) (x)- (A.22) 
Finally, for Z we obtain 

{c*°*(x),Z(y)} = {c'°\x),Z{y)}^c^O (A.23) 
Hence this Poisson bracket vanishes weakly. Considering now the Poisson bracket between c and c*°* we get 
{c'°\x),c{y)} = {c'°'{x),c{y)} + a^{x){Z,{x),c{y)} + a,{x){Z^{x),c{y)} + a{x){Z{x),c{y^ 



(A.24) 



The results of the last three individual Poisson brackets occurring above are listed below 



{Zj{x),c{y)} = {n^Jd^)q''^UaWjP){y)^^5\x,y) (A.25) 
{Z^{x),c{y)] = -{ny^d^q'''UaPS^^,){y)5\x,y) 

{Z(x),c(y)} = (^^^=]iy)6\x,y) 
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Now, we can solve equation ()A.23p for a which yields 



a(oc) = - I d V''{^k y){c'"\x),c{y)} (A.26) 



^'"^^^^ {nVd^)q'''UaPS\)^y^d^^ (x) (A.27) 



Here we reinserted the expressions for aj and derived before. The final step which includes the construc- 
tion of the Dirac bracket can again be found in the main text. 



B Comparison with Symplectic reduction 

The spatial diffeomorphism invariant quantities 

[iia),Ha)) , {f{a),P{a)) , {q,,ia),f^ (a)) (B.l) 

shown in equation (|3.3U|) are also obtained in [12j through symplectic reduction which is an alternative 
method to show that the pairs in (|3.30p are conjugate. 
To see how this works, we compute 



di 



T{a) 



d 
It 



d^x dei{dS/dx) 5{S{x),a)T{x) 



X 



J d^x det{dS/dx) (^5{S{x),a 



T{x) 



dt 



d{S{x),a) T{x) 



+ 



d r^, . 



s{xy 



+ [ d^x det{dS/dx) 



d6{a',ay 



T{x) 



a'=S(x) J 

(-5"(xR[5(5(x),c7) T(x)\ 



d 



Tlx) 



a'=S{x) 



T{x) 



S{x)=cj 

d 



d^x dei{dS/dx) 5{S{x),a) 



X 



dt 



T{x) - i-S^x)) S^{x)Ta{x) 



(B.2) 



S(x)=a 



where we have used da[Sjdet{dS/dx)] = 0. Exactly the same calculation reveals 



dt 



dt 



ax)- 

Qjk{x) 



d_ 

di^ 



S\x)] Snx)Ux) 



S{x)=cr 



d_ 

dt' 



S\x)) S^{x)qjkA^) 



(B.3) 



S{x)=cr 
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Using ()B.2p and ()B.3p we can now rewrite the symplectic potential in terms of the spatially diffeomorphism 
invariant variables as follows (.) := and J = det{dS/dx)) 



e 



X 



ilT + t P + S^ Pj+ P 



ab 



X 



ab 



(Px 



X 



ab 



d-^x 



X 



^^ + tP + S^ Pj+ q^kiSi p-') + 2q,kS^,Sl p 



J d'^x 



X 



IT 



P . S^aS'lp^'' 



J 



+ I d^xS^ 
'x 



Pj - 2da{qbcS^ p"'' 

Pj-2da{qbcS^ p" 



[^] S{.)=. + P [^] Si.)=a + fefc] 5(.)=. 



+ / d^xS^ 
'x 



d-^a 



^TT + T P + qjkP 



+ / d^x 
Ix 



Pj +SUtT e,a + PTa + p'''S'lSl,S;;qkl,a ) " 2da{qbcS^ P 



d-^a 



^ vr + T P + qjkP 



+ / d^xS^ 
Jx 



Pj + ia + PTa + p''S';iS[,da{S'kSlqef) " 25,((?b,5| p''') 



d^a 



i^ + TP + iijkP>' 



+ / d^xS^ 
Jx 



Pj + vr e,a + PTa + P''{qbc,a + 2qecSpla) " 2da{qbcS''j p"') 



+ / d^x 

IX 



d-^a 



d-^a 



i^ + fp + iijkiP'' 

Pj + (^vr e,a + PTa + -2 

l^ + fp + iijkfP'' 

in + fP + ^jkP''' 



+ / d^xS^ 

X 



QabdcP + - {2qa(b,c) - qbc,a) P 

Pj + (vr C^a + PTa + -2qabDcp'"' 



+ / d^x Sic, 

X 



,tot 
a^a 



where we used 



oa ok oe oa qH ce na Qtc ce ca oe qi 

'-'j '-',b'-'k,a — '-'j '^,ba k — '-' j '-',ab'^k — j ,b'~' ,a'-' k — 



~fa ok oe 



a ok oe oa ok oe 



(B.4) 



(B.5) 



as well as the definition of the Christoffel symbol in the second to last step (notice that p is a tensor 
density so that Dbp"'' = dbp"^ + Tl^p^""). 

Formula (IB.4[) means that on the full phase space we can switch to the new canonical pairs (IB.1[) on X 
and the canonical pair {S^ ^ Pj = 5*^0*°* = c*°*). The pairs (|B.ip are thus c*°* invariant while 5-^ is pure gauge. 
The symplectic reduction of the full phase space with respect to c*°* is therefore precisely coordinatised by 
(jB.ip which is identical to equation (j3.3Up in the main text. 
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C Effective Action and Fixed Point Equation 



The aim of the present section is to derive, at least in imphcit form, the Lagrangian that corresponds 
canonicahy to the physical Hamiltonian. This can be done by calculating the inverse Legendre transform!^ 
and leads to a fixed point equation, which can be solved order by order, in principle. Interestingly, the 
Lagrangian turns out to be local in dust time, but will be non-local in dust space. However, the Hamiltonian 
description is completely local. 

The inverse Legendre transform requires to solve for the momenta P^^{a), n(cr) in terms of the corre- 
sponding velocities Vjfc(cr), T((t), respectively, defined byEl 

Vjfc(cr) = Qjfc(cr) = {H,(5jfc((T)} 
T{a) ^ E{a) = {H,E{a)} . (C.l) 

This can be done by using the first order equations of motion for Qjkio'), S(fT), derived from the physical 
Hamiltonian H. From the physical Hamiltonian H we obtain the Lagrangian 



L[Q,F;H,T] 



dV L((j) 



dV 



-P^'V^k + ^nr] -H[Q,P;E,U] 

K \ ' 



EH 



(C.2) 



where it it is understood that the solution of (jC.ip for P-'^, H has to be inserted. 

With the dynamical lapse and shift given hy N = C/H, Nj = —Cj/H, respectively, we obtain for P^^ 



pj 



k _ \/det(Q) j^-iyjkmn^/^ _ 

2jy I [Wmn y-j^^ )mnj 



(C.3) 



with Kmn denoting the extrinsic curvature. This leads to the following expression for the velocities Vjk and 



T: 



We conclude 



Vjk = 2[N Kjk + D^jNk)] and T 



N 



Vdet(Q) 



L 
L 
L 
L 
L 
L 



dV 



-VjkP^^ + ^TH - H 



1, 
A 



(C.4) 



(C.5) 



H 



A7dit(Q) 



U-H 



N(-KikP''' + 



c 



'-K,kP^' + 



n2 



+ N'JCj-H 



^3 C 



AVdet(Q) 

n2 



Ayd^t(Q) 



C 



d^a N 



1^ \^Jde,i{Q) 



C 



dV iVVdet(Q) 



d?a N \/det(Q) 



Vdet(Q) 



\Q^%jE^k + v{E) 



+ -((V„H)^-(g^%-Hfc + KH)) 



This is possible because the Legendre transform is regular. 

Note that Qjk and H must be treated as independent variables in addition to Qjk and S. 
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In the third step we performed an integration by parts and in the fourth step we substituted the expressions 
for dynamical lapse and shift, in the sixth we rewrote P^^ in terms of Kjk, in the seventh we substituted 
for C and in the last we introduced the vector field u = jj'idr — N^d^j). 

If lapse and shift would be independent variables, the final expression in ()C.6P would coincide with the 
3+1- decomposition of the Einstein - Hilbert term minimally coupled to a Klein - Gordon field with potential 
v\ Since Nj is a constant of the physical motion and N = sjl + Q^^NjNk, we could, in particular, consider 
the case Nj = 0, whence = 1. In that case (jC.Sp would agree with the usual Lagrangian description on 
dust space-time for a static foliation. However, fundamentally lapse and shift are not independent variables, 
and we must use this fact in (jC.ip in order to solve for P^'',Il. We now turn to this task. 

By definition 



iV. 



H C H C/v/drt(Q) 
- l{D,Kf-D,Kl)+^{VuE)D,E, 



K 



Kik = — ^ = \Vik - 2D,.Nk\] , 



V,, 



^ T-Q^^NjDkE . (C.6) 



^1 + Q^^NjNk 

The set of equations (|C.6p . when inserted into each other, yields an equation of the form 

N, = G,[Nk;Qki,Vki,E,T], (C.7) 

where Gj is a local function of its arguments and their spatial derivatives up to second order (in particular, 
second spatial derivatives of Nj). Since P^^,Il are known in terms of Qjk,Vjk,E,T , once Nj (and thus A^) 
is known as a function of these arguments, we have reduced the task of performing the inverse Legendre 
transform to solving the fixed point equation (|C.7p . 

Unfortunately, (|C.7p is not algebraic in Nj, so that a solution just by quadratures is impossible. Also, 
it represents a highly nonlinear system of partial differential equations of degree two, so that linear solution 
methods fail, as well. We leave the full investigation of this system for future research. However, the fact 
that it is a system of fixed point equations suggests to look for a solution by perturbative or fixed point 
methods: 

1. If we make the Ansatz that A'^- is small, in an appropriate sense, then we may expand Gj[Nk] around 
Nk = to linear order and solve the resulting linear system of PDE's. 

2. The fixed point method suggests to write the solution in the form 

N, = G,{Gk,{Gk,{..{Gk„{..))..))). (C.8) 

If convergence is under control, then an n — th order approximation may be given in the form 

A^f ) = G,(Gfc,(G,,(..(G,„(0))..))) , (C.9) 

which consists in setting the starting point of the iteration at Nj = (which is a reasonable guess if the 
exact solution is indeed small in an appropriate sense, having a test clock in mind) and to iterate n times. 
The expression (1C.9P contains spatial derivatives of the metric of order up to 2(n + 1), but is only of first 
order in r-derivatives, thus establishing that the final Lagrangian is spatially non-local in dust space but 
temporally local in dust time. 
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D Two Routes to Second Time Derivatives of Linear Perturbations 

In this appendix we consider a general Hamiltonian system with canonical coordinates {q,p) and standard 
Poisson brackets {p, q} = 1 and a Hamiltonian function H(q,p). We will consider only one degree of freedom 
but everything generalises to an arbitrary number of degrees of freedom. 

Lemma D.l. Let (go (''"))Po (''")) be an exact solution of the Hamiltonian equations of motion 

go(r) = \^{H,q}{q,p)j^^^^^^^ 

Po{t) = [{H,p}iq,p)]^^^^^^^ (D.l) 

p=po('^) 

Define the perturbations 6q := q — qoir), 5p := p — Po{t). Let H{q,p) = Y^'^^qH^'^^ be the expansion of 
H{q,p) around qoiT),po{T) in terms of the perturbations where H^"'^ is the n-th order term in terms of the 
perturbations. Then 

(i.) Expanding the full Hamiltonian equations of motion for q, p to linear order is equivalent to using the 
function H^"^^ as a Hamiltonian for the perturbations. 

(a.) Expanding the equation for q to linear order is equivalent to the equations for 6q = 6q, 6p = 6p to linear 
order. 

Proof. Notice that qo{T),po{T) do not carry any phase space dependence in contrast to 6q,5p. Therefore 
{6q, 5q} = {5p, 5p} = and {6p, 5q} = 1. 

(i.) 

Consider the full Hamiltonian equations of motion for a general solution (9(t),p(t)) 

q{T) = HM^),p{t)) and p(r) = -i/,p(g(r),p(r)) (D.2) 

where H^q = dH/dq, H^p = dH/dp. We set dqir) = q{T) — qoir) and 6p{t) = p{t) — Po{t). Subtracting 
from ()D.2p the equations for (go(''")iPo (''")) we obtain 

6q{T) = H,p{q{T),p{T))-H,p{qo{T),po{T)) 

6p{t) = -H,,{q{T),p{T))+H,,{qo{T),po{r)) (D.3) 
which is still exact. Expanding the right hand side of (|D.3p to first order in Sq{T), Sp{t) we find 

6q{T) = H^pq{qoiT),po{T))6q{T) + H^pp{qo{T),po{T))6p{T) 

5p{t) = -H,qq{qo{T),po{T))6q{r)-H,gp{qo{T),po{r))6p{r) (D.4) 
On the other hand we have 

H^^^ = \H,,,{qo{T),po{T))[6q]^ + ^H,pp{qo{T),po{T))[Sp]^ + i/,,p(go(r),Po(r)) [5p] (D.5) 
Then it is trivial to check that ()D.4p is reproduced by 

5q{T) = , Sq} ,,^,,(^) and 6p{t) = {H^^^ , 6p} s^=8,(r) (D.6) 

(5p— (5p(t) (5p— (5p{t) 

(h.) 

Let p{t) = F{q{T), q{T)) be the solution of solving g(r) = H ,p{q{T),p{T)) for p{t). Then 

^{r) = H^pg{q{T),p{T))q{T)+H^pp{q{T),p{T))p{T) 

= H,p,{q{T),p{T))q{T) - H,pp{qiT),piT))H^,{q{r),pir)) 

= H,p,{q{T),F{q{T)A{r)))q{T) - H ,pp{q{T) , F {q{T) A{t))) H ,,{q{T) , F {qir) , q{r))) 

= : G{q{T)A{r)) (D.7) 
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Equation ()D.7p is what we mean by the q{T) form of the quations of motion, i.e an equation only involving 
q,q,q but no longer the momentum. Subtracting from ()D.7p the corresponding equation for qoir) and 
expanding the right hand side to first order we obtain with G = G{q, v) 



5q{T) = G^q{qo{T),qo{T))5q{T) + G,^(go(-r), ^o(T))(5g(r) 



(D. 



Now 



G,viq,v) 



[H,pqqiq,p)v - H^ppg{q,p)H^q{q,p) - [H^pp{q,p)H^qq{q,p)]^^p^^^^^ 

+ [H,ppqiq,p)v - H^ppp{q,p)H^q{q,p) - H^pp{q,p)H^pq{q,p)]^^p^^^^^ Fqiq,v) (D.9) 
H,pq{q,p) + [H^ppq{q,p)v - H„p{q,p)H^q{q,p) - H^pp{q,p)H^pq{q,p)]^^p^^^^^ Fv{q,v) 



Since v = H^p{q, F{q, v)) is an identity we obtain 

l = H^pp{q,F{q,v))F,{q,v) and = H^pq{q, F{q,v)) + H^pp{q, F{q,v))Fq{q,v) (D.IO) 

by taking the derivative with respect to the independent variables v, q respectively. This way we can eliminate 
the derivatives of F 



and F^q{q,v) 



H,pq{q,F{q,v)) 



Hpp{q,F{q,v)) ^ Hpp{q, F{q,v)) 

Substituting (|D.lip into (|D.9p we obtain the simplified expressions 



(D.ll) 



G,q{q,v) 



H^ppqV H pppH q 



H 



pp 



iq^p) 



p=F(q,v) 



H 



PI 



.pp 



[H,pqq{q,p)v - H^ppq{q,p)H^q{q,p) - H^pp{q,p)H^qq{q,p)]^^p^^^^^ 
-[H,ppqiq,p)v - H^ppp{q,p)H^q{q,p) - H^pp{q,p)H^pq{q,p)]^^p^^^^^ ^ 

Now we invert the second equation in ()D.4p for 5p and obtain 

. I s (^g(r) - H^pq{qo{T),po{T))5q{T) 
H,pp{qo{T),Po{T) 

Taking the time derivative of the first equation in ()D.4p and using (jP.lSP yields after some algebra 

H ^pqi^H pp H ppH pf^ 



(D.12) 



(D.13) 



H^pq H^ppH^qq 



H 



pp 



{qo{T),Po{T-)) dq{T) 



+ 



H 



,pp 



H 



pp 



{qo{T),PoiT)) 5q{T) 



(D.14) 



where e.g. H^pp{qQ{T),po{T) := ■^H^pp{qQ{T),po{T). Carrying out the remaining time derivatives in (ID.14h 
and comparing with (jP.SP evaluated with the help of ()D.12p at q = qo{T), v = qoir) = H ^p{qQ{T) , po{t) we 
see that the expressions coincide. □ 

E Constants of the Motion of n'th Order Perturbation Theory 

In this section we will show that for any fully conserved quantity F of a Hamiltonian sytem with Hamiltonian 
H, when expanding both the equations of motion and F to order n then F is still a constant of motion up 
to terms of order n + 1. 
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Let ?tio(t) = {qo{T) , po{t)) be an exact solution of a Hamiltonian system with canonical coordinates 
m = {q,p), non vanishing Poisson brackets {p,q} = 1 and Hamiltonian H = H{m) = H(q,p). Define 
6m = m — mo(T). Since mQ{T) is just a number (for fixed r) we immediately have the non vanishing Poisson 
brackets {Sp,Sq} = 1. For any function F on phase space we consider its Taylor expansion around mo{T) 
given by 



F{m) = F^''\mo{T);5m) 



(E.l) 



n=0 



where F^'^\mo{T); Sm) is a homogeneous polynomial of degree n in 6m whose coefficients depend explicitly 
on the background solution mo(r), that is 



k=0 



k\ {n-ky. 



Qbp 



[dqT [dp] 



n—k 



(mo(r)) [5qf[5p] 



n—k 



(E.2) 



Lemma E.l. The Poisson bracket {F, G} can he computed either by first expanding F,G as in iE. 1\} and 
then using the Poisson bracket for 6m or by using the Poisson bracket for m and then expanding the result 
as in IE.1\) . 



Proof. The proof is elementary: Since 

oo ^ 

F=T — 

^ k\ l\ 



k,l=0 



gk+l p 
[dqt [dpf 



(moir)) [6qt 



(E.3) 



we have with the substitution of F by 



oo 

^ kW. 



k,l=0 



Qk+l+lp 
k+1 



(mo(r)) [6qf [6pf = F,, 



(E.4) 



and similarly = F^sp- Since one computes Poisson brackets the first way by first expanding and then 
taking derivatives with respect to 6q, 6p while the second way we compute Poisson brackets with respect to 
q,p and then expand, the assertion follows. □ 

Lemma E.2. Suppose we expand the Hamiltonian to n—th order in 6m with n > 1. Suppose also that F is 

an exact constant of the motion with respect to the Hamiltonian H . Then: 

(i.) The equations of motion up to order n for 6m are generated by the Hamiltonian 



n.+l 



Hn — ^ H 



(k) 



(E.5) 



k=2 



(a.) The perturbation up to order n of F given by 



(n) 



(E.6) 



k=l 



is a constant of motion with respect to H^ up to terms of order at least n + 1. 
Notice that the Hamiltonian starts at order two and ends at order n + 1. 
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Proof. 
(i-) 

Let m{t) be any solution of the exact equation of motion. We have for example 

^(^) = [{H, q}]m=m{T) 

Subtracting the same equation for mo(r) and setting dq{T) = q{T) — qoir) we find 



(E.7) 



k=2 

from which the assertion follows immediately (the proof for 6p is identical), 
(ii.) 

Using the explicit background dependence of F^^^ = F^^\mQ{T); Smij)) we have 

dqo ddp dpo 86 q d5q ddp d6p 86 q 

(1) 



(E.J 



(E.9) 



where we used the first part of the lemma as well as the fact that H^q{mo) = H'^^^, H^p{mo) = H^l^- All 
Poisson brackets are with respect to the coordinates 6q, 6p. 
Now observe the important fact 



.90 



k ^ 
^ l\(h - 



Qk+l p 



[6q 



_d_y. 1_ 



Qk+l p 



dSq ^ (/ + 1)!((A; + 1) - (/ + 1))! [dqf+^ [8pf+'^-^'+'^ 
1 d>'+^F ...i,,,k+i-i 



'+1 rAr,i('=+i)-('+i) 



[Sqr' m 



n fc+l 

— T- 



8 
85q 



p{k+l) 



1 d'^+'F .,+1 



[6pY 



and similarly F^p^ = F^^^^\ 

Combining (|K9]) and (IKIO]) we see that 



dT 



Hence 



dF 



dT 



n \ ^ 



fc=l 



n+1 



«=2 



(E.IO) 



(E.ll) 



(E.12) 



We would like to show that the terms up to order n in (IE.12P vanish identically. Since {H^^\ F^^^} is of 
order k + I — 2, for given k we can restrict the sum over I from / = 2 until n + 2 — /c up to terms of order 
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0{5^~^^). Notice that n + 2 — A; is at least 2 (for k = n) and at most n + 1 (for k = 1) for all values of k 
which is the allowed range of /. It follows 



dK 



dT 



- +0(5"+i^ 



n+2-fc 
1=2 



E 

fc=i 

n 

E 

A:=l 
n 

E 

r=l 
r=l fc=l 



fc=l 



(E.13) 



where in the second step we have introduced the summation variable r = k + I — 2 which for given k takes 
range in k, ..,n (lowest value for / = 2 and highest value for I = n + 2 — k) whence / = r — A; + 2, in the 
third step we have changed the order of the k and r summation in the second term (keeping in mind the 
constraint 1 < k < r < n) while the summation variable k was renamed by r in the first term and in the 
fourth step we noticed that the first and second term can be combined by having the k summation extend 
to r + 1. 

Now we exploit the fact that F is an exact invariant, that is 



kJ=l 



r=0 



V+l 



^|j^(r-fc+2)^^(fe)| 



.k=l 



(E.14) 



where in the second step we collected all terms of order r (notice that l = r — k + 2> las required). 
Since ()E.13P is an identity on the entire phase space, the Taylor coefficients of [5q]^ [5pf have to vanish 
separately for all /c, / > 0. The term corresponding to order r in ()E.13p contains all terms of the form 
s = 0, ..,r. Therefore we conclude 



r+l 

Efff 

k=l 



(r-fc+2) 



(E.15) 



identically for all r. In particular, (jE.lSP implies 



dFn 
dT 



0{S 



n+l\ 



(E.16) 
□ 



The only n for which the term 0((^"~''^) vanishes is for n = 1 as one can see from ()E.12p since then 
k = 1,1 = 2 can only take one value which already contributes to order r = 1. Thus for n = 1 we even have 



dFi 
1^ 







(E.17) 



It is instructive to see how the background equations find their way into demonstrating the important result 
()E.16p which ensures that an exact invariant expanded up to order n remains an invariant up to higher 
orders for the equations of motion expanded up to order n, thus simplifying the task to integrate those 
equations of motion. 
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F Generalisation to Other Deparametrising Matter 



In this work dust was used as a reference frame in order to define a physical time evohition. We chose 
dust because then, for the case of no perturbations, the induced physical Hamiltonian yields the exact FRW 
equations used in standard cosmology. However, when perturbations of the metric and the scalar field are 
considered, deviations from the standard FRW framework occur, which are, however, still in agreement with 
observational data. Since General Relativity does not tell us which is the right clock to use for cosmology, 
we chose a clock such that the resulting physical Hamiltonian is as close as possible to the FRW-Hamiltonian 
in standard cosmology, where one uses the Hamiltonian constraint c as a true Hamiltonian. The physical 
Hamiltonian used here, i^dust = y^C^ - QvdCj, reduces to H^^^ = C for an FRW universe, namely to 
the gauge invariant version of the Hamiltonian constraint. However, the question arises how generic are the 
results obtained from a dust clock and what changes do we expect when choosing other matter than dust 
to reparametrise or even deparametrise the constraints of General Relativity. To illustrate this issue, let us 
discuss the phantom clock introduced in [9j which leads to a physical Hamiltonian of the form 

Hphan = /"f^VFphan(f^), (F.l) 

with the Hamiltonian density -ffphan defined as 



X 



where 



^phan(^T) = \l\{F{C,Ci,Q^,)) + ^-^{F{C,Ci,Qi,)f -a^Q^^{a,T)CiC,{a)Q{^^^^^ (F.2) 

F{C,Ci,Qij) := C\a,T) - Q^^(c7, r)QC7,(a) - q2Q((7, r). (F.3) 

Above we introduced the abbreviation Q = det{Qij), and a > is for the moment an arbitrary constant 
of dimension cm~'^ that enters the phantom field action as a free parameter. Recall that the expressions 
for C and Cj were given by the geometry and matter part of the total Hamiltonian and diffeomorphism 
constraint, respectively, that is 

C{a, r) = C7§^°(a, r) + C"^'^"^^(f7, r) and Cj{a) = Cf°{a) + Cf '^""''(a). (F.4) 

From now on we will drop the r dependence of C((T, r) and Qij{(J,T) in the expression for //phan and write 
them explicitly only when confusion could arise otherwise. 

For the dust Hamiltonian we saw that the first order equations of motion obtained for H, H and Qij, P^^ 
look similar to the standard cosmological equations apart from the fact that in the case of a dust clock 
we obtain a dynamical, that is phase space dependent, lapse function A'dust = C/Hdnst and shift vector 
-^dust ~ -^V-^dust- Moreover, in the general equations there occurs a term proportional to -^dust-^dusf 



With respect to Hphanj we want to analyse now whether it produces a similar effect and leads to possi- 
bly more serious deviations from the standard equations. Starting with the first order equation for H we 
obtain 



'^{ct^t) = {Hphan,S(cr,T)} 

c 




. 1 ^ {C'-Q^^aC,-a'Q){a') 

^ 4^/(1(^2 - Q'^dCj - a^QY - a^Qi^CiCjQ){a') 



(F.5) 

{C""""^'-(a'),H(a,r)} 



{C^ -Q'^CiCj+a^Q){a') 



{C~((7'),H(a,T)J 
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Introducing the dynamical lapse function A'phan and the dynamical shift covector N^^^^ as 



-^phan(o-) := 



C 



phan 



-(^) 7T + 



{C^ -QijCiCj + a^Q){a) 



phan 



(F.6) 



we can rewrite the first order equation of motion for H as 
S(a,r) = {Hphan,S(a,r)} = j dV (iVphan(aO{C-^"^^(aO, r)} + iY^hanlGr"^"^!^'), r)}) . (F.7) 



Hence, we realise that, similarly to the case of the dust clock, the effect of the phantom clock results in 
the appearance of a dynamical lapse function and a dynamical shift vector. However, due to the more 
complicated structure of Hphan as compared to Hjugt, A^phan and N^^^^ are not simply given in terms of 
C/i^phan and — Cj/i^phan respectively. Now also terms occur which include higher than linear powers of 
the constraints in the nominator and denominator. Since the quantities H and obviously Qij also Poisson 
commute with any function that does depend on Qij only, the calculation works analogously in these cases 
and we obtain 

ri(a,r) = {Hphan,n(a,r)} = J dV' (iVphan(aO{C--"^^(aO, + iV^han{C'r"'^(cT'), n(a, r)}) (F.8) 



and 

Qij{a,T) = {H^^,^,Qi,{a)} = J dV (iVphan(cTO{C'S^°(cTO, r)} + Ar^han{Cr(a'), Qi.la, r)}) . (F.9) 

X 

For the dynamical variable P*-' things look slightly different, bccaiisc P^^ does not Poisson commTitc with 
functions depending on Qij. Therefore we get, as in the dust case, an additional contribution proportional to 
the Poisson bracket {Q'^\a'), P^^{(t, t)}. Furthermore, since Hphan includes a term of the form a^Q, we also 
obtain a term proportional to the Poisson bracket {Q{a'), P^^ {a,T)}. The explicit results for these Poisson 
brackets are 



{Q''\a'),P'^{a,T)} = -{Q"'Q^^ + Q^''Q'^){a',T)6'{a',a) 
{Q{a'),P'^{a,r)} = -«(QQ^^)(a', r),5=^(a', a). 
Inserting this back into the eqn. for P^^ , we end up with 
P^ia,T) = {Hphan, P'^' (a, r)} 

= / dV (Arphan(a'){C7S^°(a'), P'^ia, r)} + N;^,,{Cr{a'), P'^{a, r)}) 



(F.IO) 



K 



2 I -f^phaniVphan-^phan 



1 

2 + 



(C2 - Qi^dCj + a^Q) 



+ 



2 



phan 



+ 



ia,T). 



(F.ll) 
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The remaining Poisson brackets in the first order equations for the dynamical variables are the same Poisson 
brackets that occur when Hjust is used as a Hamiltonian. Inserting the results obtained there into the 
corresponding equations for the case of Hpham we obtain the following final form of the first order equations: 



n(a,r) 
Qii(o-,r) 



iVphann 

Q 



(cr,r) + (£^phanH) (o-,r) 



2K 



phan 



Q 



(ct,t)(G 



where 



tjmn 



{QimQjn ~\~ QinQjm QijQ 



I] ^ mn j 1 



with its inverse given by 



(F.12) 

(F.13) 
(F.14) 



For the gravitational momentum we have 

Qmi 



phan 



^/ditQ 

+ V^ditQ[G-i]^'=™"( (D^D^iVpi: 



-^phan 

n -1 



iVphaniim„[Q] ) + ^ iVphan V^ditQ H'^H'^' + {C^^^^Py'' ) t) 



2 



+- 



-f^phan^^'phan-^phan 



1 



+ 



phan 



+ 



(a,r). 



(F.15) 



One can see that the first order equations for H, H and Qij in equation ()F.12p are identical to those for 
the dust Hamiltonian Hdust, apart from the different definitions of the dynamical lapse function A^phan 
and shift vector iV^^an ™ equation ()F.6p . However, the equation for P*-' differs from the corresponding 
dust Hamiltonian equation. The term in the second last line in equation (jF.lSp corresponds to the term 
— |(i7dust-^dust-^dust)('^) ™ equation for P^^ derived from H^ust- In the present case this term looks a 
bit more complicated, since we have to divide the whole expression by the term in the square brackets which 
is identical to one in case of Hjusf The additional term in the last line of equation (jF.lSP comes from the 
terms a^Q and a^Q"^^ CiCj in Hphan, which are absent in Hjust- Since in this term we cannot factor out C or 
Ci but only , we are also not able to reexpress this term by means of the dynamical lapse function A'^phan 
and the dynamical shift vector A^p^an' respectively, as it was possible for the term in the second last line. 

In summary, when using a phantom scalar field as a clock we also obtain deviations from the standard 
treatment in which the Hamiltonian constraint is used as a true Hamiltonian. These deviations manifest 
themselves in the appearance of a dynamical lapse function A^phan and a dynamical shift vector A^p^an' 
ogous to the case where dust is used as a clock. However, the explicit dependence of A^phan and A^p^an ^^e 
dynamical variables is more complicated than for the corresponding quantities N^ust and N^^^^. Another 
modification occurs in the equation for the gravitational momentum P^^ . While it contains a term that 
is second order in A^p^an' ™ complete analogy with the case of H^ust) it also features an additional term 
proportional to a^. 
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The Special Case of an FRW Universe 

It is interesting to study the special case of FRW also for Hphan- Recall from [15] that by assuming 
homogeneity and isotropy, Hdustj-^dust and A^'^^s^ reduce to the following quantities 

iVdust(^)^l, K,,,ia)'-^0, i7d,.t(^) ^ C--(a), (F.16) 



where 

CirW = - 6(^)' + 2A) + ^ (E^ + V{E)) )(r). (F.17) 



Here A = A(t) is a function of dust time r and the dot refers to a derivative with respect to dust time. In 
particular A{t) can be understood as the gauge invariant extension of the ordinary scale factor a{t) used 
in standard cosmology. The difference between the two is that A is gauge invariant and thus a physical 
observable whereas a is not, since it does not commute with the Hamiltonian constraint of FRW. 
The gravitational canonical variables are given by Qij = A'^{T)5ij and P^^ = —2A{t)6^^. We mentioned 
previously that a consequence of this behaviour is that the unperturbed equations of motion for (H, H) and 
{Qij, P^^) agree with the FRW equations used in standard cosmology. In particular, the deviation from the 
general standard equation of motion for P^^ vanishes in the case of FRW, because it is quadratic in N^^^^. 
For the phantom Hamiltonian Hphan things look slightly different. Here we have the following behaviour of 
-f^phari) -^phan and N^i^^^, when a homogenous and isotropic universe is considered: 



FRW 



^ -"phan ' ^ 



(F.18) 



with 



C;Cnir) = J^;^{r){l{-G{j)\2Ai^^^^^^ (F.19) 

Hence, the dust clock and the phantom scalar field clock agree only if the parameter a is chosen to be tiny 
compared to the Hamiltonian constraint Cp^^, e.g. aA^ <C Cp^^. Consequently, the equations of motion 
generated by Hphan also deviate from the standard FRW equations. The significance of this deviation again 
depends on the specific value of the parameter a, as was discussed in detail in [9]. For completeness we also 
list the first order equations of motion here 



■(r) = {N;Cn-^){r) (F 



n 

n(r) = -{N;Cn^V'{E))ir) 
2 

Qij = {^ph^-^^GijmnP"^"){T) 



P'' = KCnir) {2P^rnprn _ p^, pmn^ ^ ^Q^^ C^Z. " Q^^VQ(2A + ^V{E)) + ^^Q*^) 

Taking into account that Qij = A^5ij and using the first order equation for Qij, we solve for the momenta 
P*-^ = — 2^/A'j^]^j^(5*-' in terms of Qjk = 2AA. In order to derive the corresponding FRW equation with 

respect to the time generated by Hp^^, we take the time derivative of the equation for Qij and insert into 
the resulting equation for Qij the expression for P^^ and P^^ given above. This yields 

(i) = + - rx^¥ - ^ wrn)V(H)) + - ^ wrn)^(^). (F.21) 
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Apart from the lapse functions in the equation above which are not present in the standard FRW equations, 
we get an additional term including the time derivative of the lapse function. Using the explicit definition 
of the lapse function, we can perform this time derivative, leading to 



A^FRw \aI ('/v™w\2 \ a) ' (F.22) 

^^phan W^phan^ 

Consequently, equation ()F.2ip can be rewritten as 

The term (A/A)'^ can be replaced by considering the energy conservation law H^^^ = , that is H^^^^ = eo, 
from which we get Cp^^ = + o?A^. Solving this equation for (A/yl)^ yields 



= (A^pTrn)'A + ^(iVp\Tn)' (Pmatter + Pphan) , (F.24) 

where we used in the last line 



Pmatter = \ ^^l^ + ^^(^) ^^^^ Pphan = 1 + (F.25) 

Reinserting equation ()F.24p into equation ()F.23p . we obtain the phantom FRW equation given by 



(F.26) 

whereby we introduced 

Pmatter = 2 (A^™^)2" ~ 2^^'^'' 

1 d / 3 \ O? A? 

Pmatter = ~ Y^~dA\ ^^"^^ ~ 'C^^ ' (F.27) 

That this equation agrees with the one derived in [9j can be seen when expressing {J^^^^)^ in terms of the 
deviation parameter x := e^jc? A^ used there, resulting in (-/Vp^^)^ = 1 + l/x. 

In general, choosing one clock or the other might have significant effects on the equation of motion. Gen- 
eral Relativity does not tell us which clock is convenient to work with, hence additional physical input is 
needed. The results of the application of this framework for FRW in [T5] show that choosing dust as clock 
reproduces the standard FRW equations. Thus we could call the dust clock the FRW-clock. Since so far 
an (approximate) FRW universe is in agreement with observational data, dust seems to be a good choice. 
However, the a parameter in Hp[J^ can be chosen such that the resulting equation of motion also do not 
contradict present experiments. Therefore, based on experimental constraints, none of the two clocks is 
excluded, nor is one of them preferred. From a theoretical point of view, the choice of a clock is mainly 
guided by the requirement that the constraints can be deparametrised, that is, they can be written in the 
form C = p'^io'^k _j_ ^ciock_ jjgj-g Hamiltonian density //^lock t^^^^ depend on the clock variables 
anymore, and furthermore it should be positive definite. Additionally the structure of H^'^'^^^ should not be 
too complicated such that calculations of, for instance, the equation of motions are still possible. 
However, in principle, we have a large amount of freedom to choose a clock, as long as the induced equations 
of motion do not contradict experiments. 



1 \ / \ /I 
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G Linear Perturbation Theory: Some Calculations in More Detail 



In section 16.21 we derived the second order equation of motion for the Hnear perturbation of the (manifestly) 
gauge invariant three metric SQjk- For that we needed the perturbation of the geometry and matter part 
of the gauge invariant Hamiltonian constraint, denoted by C^'^° and C"^^^^^^ , respectively, as these terms 
occur in the third term on the right-hand side of the unperturbed equation of motion for Qjk, equation 
(j4.24p . We omitted the details in the main text due to their length, and also because it turns out that 
several terms cancel when inserted back into the expression of the perturbation of the third term, shown in 
equation (j6.20p . For the interested reader, however, the detailed perturbations of the constraints are given 
below. 

The perturbed geometry constraint 6C^'^" is given by 
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Here we used that the perturbation of the Ricci scalar can be written as 
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For the perturbed matter part of the constraint SC^^^^'^^ we obtain 

^^matter 



1 /7=r7fc A'^^"' A^'^ \ ^matter I N'' /-^matter \/ dct Q f-^jk— — 



N 



+ 



\/ det Q—j 



■jm-p=rkn I 



Q' Q 



+ 



Vdet Q 

2 ' 

NX 



■jk- 



y/detQ / 1 



V_detQ— jm-^jfc^ ^ 

2^ v ^,m^,n 



6Q 



jk 



N 



■^matter \/ det Q /-^jk— — , 

C — \Q ^jr.^k + v[ 



2X 



d_ 



N 



6E. 



(G.3) 



6N, 



74 



Since the perturbation of 6C = 5C^'^° + 5(7™^**'''' occurs in equation (j6.20p multiplied by a factor of _ , 
we will present it here already with this factor in front 
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Going back to the second order equation of motion for Qjk shown in equation ()4.24p . we remind the reader 
that the perturbation of the first term on the right-hand side involves a term that we had already calculated 
for the equation of motion of 5E. For this reason, we presented in the main text only the perturbation of 
the remaining term {Qjk — {C^Q)jk), not the final result for the full first term. For those interested in more 
detail, we display it here: 
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H Gauge Invariant Versus Gauge Fixed Formalism 



In this appendix we investigate the question under which circumstances a manifestly gauge invariant for- 
mulation of a constrained system can be equivalently described by a gauge fixed version. 

We begin quite generally and consider a (finite dimensional) constrained Hamiltonian system with first 
class constraints Cj, I = l,..,m on a phase space with canonical pairs {q"',Pa), o, = m < n. If 

there is a true, gauge invariant Hamiltonian H (not constrained to vanish) enlarge the phase space by an 
additional canonical pair {q^ ,po) and additional first class constraint Cq = po + H. The reduced phase space 
and dynamics of the enlarged system is equivalent to the original one, hence we consider without loss of 
generality a system with no true Hamiltonian. 

The canonical Hamiltonian of the system is a linear combination of constraints 

^can = X^Cl (H.l) 

for some Lagrange multipliers whose range specifies the amount of gauge freedom. A gauge fixing is defined 
by a set of gauge fixing functions Gj with the property that the matrix with entries M/j := {Ci,Gj} has 
everywhere (on phase space) non vanishing determinant!^. Notice that we allow for gauge fixing conditions 
that display an explicit time dependence. The conservation in time of the gauge fixing conditions 

= = ^Gl + {Fean, Gl} = ^Gl + \^ Mji (H.2) 

at at at 
uniquely fixes the Lagrange multipliers to be the following phase space dependent functions 

A^ = -^(M-y^=:AS (H.3) 

By arbitrarily splitting the set of canonical pairs {q"',Pa) into two sets (T^ , vr/ ) , / = 1 , . . , m and {Q^,Pa), A = 
1, ..,n — m we can solve Cj = Gj = for 

C'j = TTi + h'j{Q,P)=0, G'j = T^ -T^{Q,P) = (H.4) 

for certian functions h, r which generically will be explicitly time dependent. The variables T, vr are called 
the gauge degrees of freedom and Q, P are called the true degrees of freedom (although typically neither of 
them is gauge invariant). 

The reduced Hamiltonian H^ediQ,P), if it exists, is supposed to generate the same equations of motion 
for Q,P as the canonical Hamiltonian does, when the constraints and the gauge fixing conditions are satisfied 
and the Lagrange multipliers assume their fixed values ()H.3p . that is, 

{-f^red, /} = {-f^can, f}c=G=X~Xo=0 = [^o{Cl, /}]c=G=A-Ao=0 (H.5) 

for any function / = f{Q,P). For general gauge fixing functions the reduced Hamiltonian will not exist, 
the system of PDE's to which (jH.SP is equivalent to, will not be integrable. 

However, a so called coordinate gauge fixing condition G^ = — with independent of the phase 
space always leads to a reduced Hamiltonian as follows: We can always (locally) write the constraints in the 
form (at least weakly) 

Ci = Mij{7:i + hi{T,Q,P)) (H.6) 
Then, noticing that M/j = {Ci,Tj}, (flL5l) becomes 

{Fred, /} = KMuihi, /}]c=G=A-Ao=o = /}]g=o = /} (H.7) 



'Ideally, the gauge Gj = should define a unique point in each gauge orbit. 
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with hj = hj{T = T, Q, P) and we used that / only depends on Q, P. This displays the reduced Hamiltonian 
as 

H,M, P; t) = fi{t)hj{T = r(t), Q, P)) (H.8) 

It will be explicitly time dependent unless fj is time independent and hj is independent of T, that is, unless 
those constraints can be deparametrised for which f/ ^ 0. Hence, deparametrisation is crucial for having a 
conserved Hamiltonian system. 

On the other hand, let us consider the gauge invariant point of view. The observables associated with 
f{Q, P) are given by 

Of{T) = [expiP'Xi) ■ /]/j=.-T (H.9) 

where we have denoted the Abelian Hamiltonian vector fields Xf by Xj := {ttj + hj, .}. Consider a one 
parameter family of flows 1 1— > T^(t) then with Of{t) := Of{T{t)) we find 

w J^L rJi RJn 

^Ofit) = f\i) XiXj,..Xj^^ . f (H.IO) 

n=0 

On the other hand, consider Hj{t) := Ohj{T{t)), then 

{Hj{t),Of{t)} = =0{,,^jyiT{t))=0{^^jy{T{t)) = Oxrf{r{t)) 

= ^ XiXji-Xjn-f (H.ll) 

n=0 

where in the second step we used that neither hj nor / depend on vrj, in the third we used that / does not 
depend on T"^ and in the last we used the commutativity of the Xj. Thus the physical Hamiltonian that 
drives the time evolution of the observables is simply given by 

H{t) := f'{t)hj{T{t), Ogit), Op{t)) (H.12) 

This is exactly the same as (IH.SP under the identification / <^ Oj(0). Hence we have shown that for suitable 
gauge fixings the reduced and the gauge invariant frameworks are equivalent. Notice that it was crucial in 
the derivation that {T^,ttj) and {Q^,Pa) are two sets of canonical pairs. If that would not be the case, 
then it would be unclear whether the time evolution of the observables has a canonical generator. 

The power of a manifestly gauge invariant framework lies therefore not in the gauge invariance itself. 
Rather, it relies on whether the gauge fixing can be achieved globally, whether it can be phrased in terms 
of separate canonical pairs, whether the observer clocks Tj are such that reduced Hamiltonian system is 
conserved and whether they do display the time evolution of observables as viewed by a realistic observer. 
In particular, the reduced Hamiltonian by construction only depends on {Q^, Pa)- In application to gravity, 
if one would not add an additional matter component such as the dust, then the true degrees of freedom 
would sit in two of the six canonical pairs {qabiP"'^)- However, notice that it is not obvious how to split 
these pairs into gauge and true degrees of freedom in an at least spatially covariant way and moreover it is 
not possible to solve for four of the p"'^ algebraically because the spatial diffeomorphism constraint involves 
their derivatives. Hence the physical or reduced Hamiltonian would become non-local. Furthermore, if one 
uses gravitational degrees of freedom for reduction then it is clear that onbe does not get the full set of 
Einstein's equations as evolution equations which is something that one may want to keep. Finally, the 
reduced Hamiltonian will not reduce to the standard model Hamiltonian in the flat space limit (i.e. with 
unit lapse) nor will it be necessarily positive. Of course, when adding matter like our dust, then similar 
to the Higgs mechanism the four dust degrees of freedom get absorbed by the metric which develops four 
additional Goldstone modes. These modes should decouple and they do as we showed explicitly in this paper 
because of the existence of an infinite number of conserved charges, however, it is not granted to happen 
when adding arbitrary matter. 
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We close this section by verifying that the reduced Hamiltonian for the Brown - Kuchaf dust model with 
the obvious choice for the gauge degrees of freedom indeed agrees with the physical Hamiltonian. As gauge 
fixing conditions we choose 

G{x) = T{x) - T{x;t), G^{x) = S^{x) - (H.13) 

whence {t, x) = (x) is not explicitly time dependent. The stability of ()H.13P with respect to the canonical 
Hamiltonian 

Hc.n= I cfx{n[c- ^/P^T^f^b]+n''[PTa + PjSi + Ca]} (H.14) 

Jx 

fixes lapse and shift to be 



no = , no = (H.15) 

Hence for any function / independent of the dust degrees of freedom 

{^^cani /}c'°'=c'°'=n-no=n-no=0 
= / d^X (no [{C, /} , ^ ^ ^ {q'^'^CaCb, /}])ct°t=0 

= 1^ d'xf^[{c,f}-^^{q'''caCbJ}] 

= / d\f{hj}] (H.16) 

JX 

where we used 



c*°* = ^ -P = h= V^^^q^^CaCb (H.17) 
Thus the reduced Hamiltonian for f = 1 equals the physical Hamiltonian under the identification qab = 



References 

[1] H. Noh and J. Hwang, Second-order perturbations of the Friedmann world model, Phys. Rev. D69, 
104011 (2004). 

[2] K. Nakamura, Gauge invariant variables in two-parameter nonlinear perturbations, Progr. Theor. Phys. 
110, 723 (2003). 

[3] K. Nakamura, Second-order gauge invariant perturbation theory, Progr. Theor. Phys. 113, 723 (2005). 
[4] K. Nakamura, Second-order gauge invariant cosmological perturbation theory, Progr. Theor. Phys. 117, 

723 (2007). 

[5] C. Rovelli. Loop quantum gravity. Living Rev. Rel. 1 (1998), 1. [gr-qc/9710008 

T. Thiemann. Lectures on loop quantum gravity. Led. Notes Phys. 631 (2003), 41-135. gr-qc/0210094] 

A. Ashtekar and J. Lewandowski. Background independent quantum gravity: a status report. Class. 

Quant. Grav. 21 (2004), R53. |;gr-qc/0404018 1 
[6] P. A. M. Dirac. Lectures on Quantum Mechanics^ (Belfer Graduate School of Science, Yeshiva University 

Press, New York, 1964). 

[7] V. Mukhanov, H. Feldman and R. Brandenberger. Theory of cosmological perturbations. Part 1. Clas- 
sical perturbations. Part 2. Quantum theory of perturbations. Part 3. Extensions. Phys. Kept. 215 
(1992) 203-333. 

[8] Conceptual problems of quantum gravity. A. Ashtekar and J. Stachel (eds.). Based on proceedings of 
the 1988 Osgood Hill Conference, (Boston, Basel, Berlin, Birkhaeuser, 1991) 

[9] T. Thiemann. Solving the problem of time in general relativity and cosmology with phantoms and 
K-essence. [astro-ph /0607380] 



78 



[10] C. Rovelli. What is observable in classical and quantum gravity? Class. Quantum Grav. 8 (1991), 
297-316. 

C. Rovelli. Quantum reference systems. Class. Quantum Crav. 8 (1991), 317-332. 

C. Rovelli. Time in quantum gravity: physics beyond the Schrodinger regime. Phys. Rev. D43 (1991), 

442-456. 

C. Rovelli. Quantum mechanics without time: a model. Phys. Rev. D42 (1990), 2638-2646. 
[11] B. Dittrich. Partial and complete observables for Hamiltonian constrained systems. Gen. Rel. Grav. 39 
(2007) 1891 (gr-qc/0411013i . 

B. Dittrich. Partial and complete observables for canonical general relativity. Class. Quant. Crav. 23 

(2006),6155-6184. [gr-qc/0507106 1 
[12] J. Brown and K. Kuchaf. Dust as a standard of space and time in canonical quantum gravity. Phys. 

Rev. D51 (1995), 5600-5629. [ gr-qc/9409001| 
[13] D. Langlois. Hamiltonian formalism and gauge invariance for linear perturbations in inflation. Class. 

Quant. Crav. (1994) 11 389-407. 
[14] M. Iftime and J. Stachel. The hole argument for covariant theories. Cen. Rel. Crav. 38 (2006) 1241-1252. 
[15] K. Giesel, S. Hofmann, T. Thiemann and O. Winkler. Manifestly gauge-invariant general relativistic 

perturbation theory: II. FRW Background and First Order. [arXiv:0711.0117 [gr-qc]] 
[16] N. Straumann. From primordial quantum fluctuations to the anisotropics of the cosmic microwave 

background radiation. Annalen Phys. 15 (2006) 701-847. |hep-ph/0505249| 
[17] V. Mukhanov. Physical foundations of cosmology, (Cambridge, Cambridge University Press, 2006). 
[18] V. I. Arnold. Mathematical methods of classical mechanics, (New York, Berlin, Heidelberg, Springer- 

Verlag, 1989). 

[19] K. Giesel, S. Hofmann, T. Thiemann and O. Winkler. Manifestly gauge-invariant general relativistic 
perturbation theory: III. FRW background and Second Order. In preparation. 

[20] T. Thiemann. Reduced phase space quantization and Dirac observables. Class. Quant. Crav. 23 (2006), 
1163-1180. [gr-qc/0411031j 

[21] C. Rovelli. Quantum Gravity, (Cambridge University Press, Cambridge, 2004). 

[22] Thiemann. Modern Canonical Quantum General Relativity, (Cambridge University Press, Cambridge, 
2007). [gr-qc/0110034| 

[23] K. Giesel and T. Thiemann. Algebraic quantum gravity (AQG) I. Conceptual setup. Class. Quant. 
Crav. 24 (2007) 2465-2498. (gr-qc/ 0607099 1 

K. Giesel and T. Thiemann. Algebraic quantum gravity (AQG) II. Semiclassical analysis. Class. Quant. 
Grav. 24 (2007) 2499-2564. [gr-qc/0607100] 

K. Giesel and T. Thiemann. Algebraic quantum gravity (AQG) III. Semiclassical perturbation theory. 

Class. Quant. Grav. 24 (2007) 2565-2588. igr-qc/0607101 1 
[24] K. Giesel and T. Thiemann. Algebraic quantum gravity (AQG) IV. Reduced phase space quantisation 

of Loop Quantum Gravity. [arXiy:0711.0119 [gr-qc]] 
[25] B. Dittrich and J. Tambornino. A perturbative approach to Dirac observables and their spacetime al- 



gebra. Class. Quant. Grav. 24 (2007) 757. arXiv:gr-qc/0610060 



B. Dittrich and J. Tambornino. Gauge invariant perturbations around symmetry reduced sec- 
tors of general relativity: Applications to cosmology. Class. Quant. Grav. 24 (2007) 4543-4586. 
[ arXiv :gr-qc/0702093). 

[26] R. M. Wald. General Relativity, (The University of Chicago Press, Chicago, 1989). 

[27] R.K. Sachs, in Relativity, Groups and Topology, ed. B. De Witt & C. De Witt (Gordon&Breach, London, 
1964), p. 523. 

[28] J.M. Stewart and M. Walker, Perturbations of spacetime in general relativity, Proc.R.Soc. London A341, 
49 (1974). 

[29] M. Bruni, S. Matarrese, S. Mollerach and S. Sonego, Perturbations of spacetime: gauge transformations 
and gauge invariance at second order and beyond. Class. Quantum Grav. 14, 2585 (1997). 

[30] S. Sonego and M. Bruni, Gauge dependence in the theory of non-linear spacetime perturbations, Com- 
mun. Math. Phys. 193, 209 (1998). 



79 



[31] M. Bruni and S. Sonego, Obscrvablcs and gauge invariance in the theory of nonhnear spacetime per- 
turbations, Class. Quantum Grav. 16, L29 (1999). 

[32] G. F. R. Elhs and M. Bruni, Covariant and gauge-invariant approach to cosmological density fluctua- 
tions, Phys. Rev. D 40, 1804 (1989). 

[33] M. Bruni, P. K. S. Dunsby and G. F. R. Ellis, Cosmological perturbations and the physical meaning of 
gauge-invariant variables, Astrophys. J. 395, 34 (1992). 

[34] C. Clarkson, Cosmological density fluctuations and gravity waves: A covariant approach to gauge- 
invariant nonlinear perturbation theory, Phys. Rev. D 70, 103524 (2004). 

[35] D. Langlois and F. Vernizzi, Evolution of nonlinear cosmological perturbations, Phys. Rev. Lett. 95, 
091303 (2005). 

[36] M. Henneaux and C. Teitelboim. Quantization of Gauge Systems, (Princeton University Press, Prince- 
ton, 1992). 

[37] R. Beig and B. G. Schmidt. Einstein's equations near spatial infinity. Commun. Math. Phys. 87 (1982), 
65. 

[38] N. Bartolo, S. Matarrcsc and A. Riotto, The full second-order radiation transfer function for large-scale 
CMB anisotropics, JCAP 05, 010 (2006). 



80 



